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Abstract 

A general pedagogical survey of the basics of the Bethe-Salpeter equation and in 
particular Wick-Cutkosky model is given in great technical details. 
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Bethe-Salpeter amplitude 

The Bethe-Salpeter approach to the relativistic two-body bound state prob- 
00 ■ lem assumes that all needed information about a bound state \B > is con- 

pu,: 



tained in the B-S amplitude 

$(x h x 2 ;P B ) =< 0\T(f> 1 ( Xl ) ( f )2 (x 2 )\B > (1) 
and in its conjugate 



><: ^{x ll x 2 -,P B )=<B\T(f>i{x 1 )<f>i{x2)\0> ■ (2) 

Here Pg is the bound state 4-momentum and field operators are in the Heisen- 
berg representation. 
Note, that 

< B\T<f>t(xi)<l>i(x2)\0 >=< OlffaMfafa^B >* 

T being antichronological operator. So we can say that $ is obtained from 
$ through time reversal. 

It is clear that instead of individual Xi,x 2 coordinates more useful are 
variables which characterize two particle relative motion and the motion of 
the system whole. 

The coordinate corresponding to the relative motion is of course x = 
X\ — x 2 . Then any linear independent combination X = T]\X\ + 772^2 can 
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serve for a description of the system as a whole. Linear independence means 
iji + / 0. In order to have a correspondence to the nonrelativistic center 
of mass definition, we take 

mi m 2 

m = — , — , m = — ; • 3) 

TYl\ + m,2 TTli + 777-2 

X\ = X + r]2X , X2 = X — T}\x equalities and the fact that P 4-momentum 
operator is the space-time translation operator enables us to write 

Mxi) = e iPx Mmx)e~ iPx , 2 (x 2 ) = e iPx <h(-mx)e- iPx • 
Inserting this expression into (1) we get 

< O|T0i(xi)0 2 (x 2 )|5 >= O(x ) < O|0i(xi)0 2 (x 2 )|B > + 
6(-x o )<O|0 2 (x 2 )0i(xi)|y3>= 
S(xo) < O|e^0i(77 2 x)0 2 (-77ix)e-^|y3 > + 
O(-xq) < 0\e iPx (f)2(-mx)(f>i(mx)e- iPx \B >= 
e- lPBX {G(x ) < O|0i(77 2 x)0 2 (-77ix)|y3 > + 
G(-x ) < O|0 2 (-77ix)0i(?72x)|B >} = 
e -iP B x < o\T ( f >1 {r }2 x)(h{-r}ix)\B > . 
So, if the reduced B-S amplitude is introduced 

p B ) = (2tt) 3 / 2 < o\TMm^M-m^\B > (4) 

when 

$(x u x 2 ] P B ) = {27r)-V 2 e- lPBX ${x; P B ) . (5) 

Analogously 

$(x b x 2 ; P B ) = (27r)- 3 /¥ p * x <i(x; P B ) , (6) 

where 

P fl ) = (2tt) 3 / 2 < y3|r0+(77 2 x)0+(-77ix)|O > . 

An equation for the B-S amplitude can be obtained with the help of 4-point 
Green's function. 
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Equation for Green's function 



Let us consider 4-point Green's function 

G(x h x 2 ; yu V2) =< O|T0i(xi)0 2 (^2)0i"(?/i)0f(?/2)|O > . 

Field operators here are in the Heisenberg representation, so we have the 
full Green's function, which can be expanded into a infinite perturbation 
series. Let us rearrange the terms of this series: first of all we sum up self- 
energy insertions in propagators of 0i and 02 particles, this will give us the 
full propagators for them, then from the remaining diagrams we separate 
(01 + 4>2) -two-particle irreducible ones, the sum of which will play the role 
of interaction operator (a diagram is (0i + 02 ) -two-particle reducible, if by 
cutting one 0i and one 02 of inner lines it can be divided into two disconnected 
parts). More vividly this is illustrated graphically in Fig.l. 




Figure 1: Integral equation for the full Green's function. 



Here / is the sum of (0i + 02)-two-particle irreducible diagrams without 
external propagators, as is shown in Fig. 2. 




Figure 2: Integral equation kernel. 



The analytical expression of the above integral equation (Fig.l) looks like 

G{x h x 2 \ yi, y 2 ) = Ai(xi - y 1 )A 2 (x 2 - y 2 ) + 
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J dzi dz 2 dz[ dz' 2 /\i{xi - zi)A 2 (x 2 - z 2 )I(zi, z 2 ] z[, z 2 )G(z[, z 2 ; y h y 2 ) . (7) 

We get the so called ladder approximation if the full propagator A(x — y) 
is replaced by the free propagator and in the interaction function / only the 
first single-particle-exchange term is left: 




Figure 3: Ladder approximation. 



An appearance of the iterative solution of this equation explains the origin 
of the approximation name: 




Figure 4: Iterative solution in the ladder approximation. 



Momentum representation 

Eq.7 is written in the coordinate space, as is well known, a relativistic particle 
can't be localized in the space-time with arbitrary precision. Therefore in 
relativistic theory momentum space is often more useful than coordinate 
space. 

In order to rewrite Eq.7 in the momentum space, let us perform the Fourier 
transformation 

G(x h x 2 ;y h y 2 ) = (2tt)- 8 / dpdqdP e~ ipx e i(iy e- iP ^- Y ^G{p,q-P) . (8) 
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Because of the translational invariance, G depends only on coordinate differ- 
ences. In the role of independent differences drawn up from the xi, x 2 , yi, 1/2 
coordinates, we have taken x = x\ — x 2 , y = y\ — y 2 and X — Y = (ijixi + 

let us clear up some technical details of the transition to the momentum 
space. First of all Ai(x\ — yi)A 2 (x 2 — y 2 ) product should be rewritten in the 
form of Eq.8. Momentum space propagator is defined through the Fourier 
transform 

AO) = (2tt)- 4 / dpe~ ipx A(p) . 

Then 

A1O1 - yi)A 2 {x 2 - y 2 ) = 

(2tt)- 8 J dp 1 dp 2 e- ipi ^- y ^e- ip ^- y ^A 1 {p 1 )A 2 {p 2 ) . 

Let us make pi = T}\P + p 2 = i} 2 P — p variable change in the integral with 
D ^ P p P p^ =771 + 772 = 1 unit Jacobian 

Ai(xi - y 1 )A 2 (x 2 - y 2 ) = 

(27r)- 8 / dPdpe- lP{x - Y) e- ipx e ipy A 1 (r h P+p)A 2 (r ]2 P-p) , 
which, it is clear, can be rewritten also as 

A1O1 - yi)A 2 (x 2 - y 2 ) = 

(2tt)- 8 J dP dpdqe- iP{x - Y) e- ipx e iqy 5{p - q)A 1 (r) 1 P + p)A 2 (r] 2 P - p) . 
Now let us transform A\A 2 IG integral 

A X A 2 IG = 

J dzidz 2 dz[dz' 2 Ai{xi - zi)A 2 {x 2 - z 2 )I{z x , z 2 ] z[, z' 2 )G(z' 1 , z' 2 \ y 1 , y 2 ) = 
(2tt) ~ 24 J dz x dz 2 dz' x dz' 2 dP dp dq dP' dp' dq' dP" dp" dq" x 
e -iP(x-z) e -i P x e iq>z § ( p _ q^A^P + p)A 2 (r] 2 P - p) x 
e- iP ^ z - z ^e- ip ' z e^ z 'l(p\ q"; P')e~ iP "~ Y ^ ] e~ ip " z ' >G(p", q; P") = 
(2tt) - 2A J dz dz' dZ dZ' dP dp dq dP' dp' dP" dp" dq" x 

AifaiP + p)A 2 ( m P - p)I(p', q"; P') x 
G(p'\ q- p»y-^-v) e -iz\v"-4<) e -iZ[P<-P) e -iZ\P"-n . 
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Coordinate integrations will produce ^-functions and some 4-momentum in- 
tegrals become trivial. At the end we get 

A 1 A 2 /G = 

(2tt)- 8 / dPdpdqdq , A 1 {r h P + p)A 2 {r ]2 P -p)I(p,q'; P)G{q',q; P) . 
After all of these, Eq.7 easily is rewritten in the momentum space 

G(p, q; P) = 5(p - q)A 1 (r h P + p)A 2 { m P - p)+ 
A 1 (r ]l P+p)A 2 (r ]2 P-p)Jdq , I(p,q , ;P)G(q , ,q;P) , 

or 

[AifeP + p) A 2 ( m P - p)]- l G{p, q; P) = 

5(p -q)+j dq' /(p, q'- P)G(q' \ q- P) . (9) 

If the following definitions are introduced 

(A ■ B)(p, q;P) = J dq' A(p, q'; P)B(qf, q; P) , 

K(p,q;P) = [A 1 (r h P + p)A 2 (r ]2 P-p)}- 1 8(p-q) , 
this equation takes the form 

K>G = 1 + I>G. (10) 

Two-particle bound states contribution 
into Green's function 

Inserting £ |n >< n\ = 1 full set of states, the Green's function can be 
rewritten as 

G(x h x 2 ;y h y 2 ) =< O|T0i(xi)0 2 (x 2 )^(?/i)0f(?/2)|O >= 

e(min[(xi) , (ar 2 )o] - max[(yi) , (2/2)0]) x 

< 0|T<fc(zi)&(z2)|n >< ^|T0+( 2/1 )0+( 2 / 2 )|O > + 

e(min[(xi) , (1/1)0] - max[(x 2 ) , (2/2)0]) x 

< 0\T<t>!(xi)<t>t(yi)\n >< n|T0 2 (x 2 )0+(2/ 2 )|O > + 

e(min[(xi) , (1/2)0] - max[(x 2 ) , (2/1)0]) x 

< 0\T<f> 1 (x 1 )<f4(y 2 )\n >< n|T0 2 (x 2 )^(2/i)|O > + 
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e(min[(x 2 )o, (2/1)0] - max[(rci) , (2/2)0]) x 

< O|T0 2 (x 2 )0+( 2/l )|n >< n|T0i(a;i)^-(^)|O > + 

e(min[(x 2 )o, (2/2)0] - max[(xi) , (2/1)0]) x 

< O|T0 2 (x 2 )0+(2/ 2 )|n >< n|T0 1 (xi)0+(2/i)|O > + 

9(min[(2/i) , (2/2)0] - max[(xi) , (x 2 ) ])x 

< O|T0+(2/i)0 2 f (2/2)|n X nlT^^O^^IO > (11) 

To separate the contribution from two-particle bound states \B >, let us 
substitute 

£|n >< n\ — > J dP5{P 2 - m 2 B )O{P )\B >< B\ . 

(< P\Q >= 2Pq5(P — Q) normalization is assumed). 

An action of two annihilation operators over \B > is needed to obtain 
the vacuum quantum numbers. Only in this case we get a nonzero matrix 
element in (11). Therefore two-particle bound state contribution into the 
Green's function is 

B = JdP < Q\Tcj )1 {x 1 )(j) 2 {x 2 )\B >< B\T<f>t (2/1)^ (2/2) |0 > x 

5(P 2 - m|)e(P o )0(min[(xi)o, (x 2 )o] - max[(2/i) , (2/2)0]) • 
dPo integration can be performed by using 

5(P 2 - m|)B(P ) = 

e(p ) 



2P 







S(P - \jP 2 + m 2 B ) + 8(P + V^ 2 + ™l)j = 

e(p ) 



2P 

and remembering (5) and (6) : 



S(P - V p2 + m%) 



B = (2tt)- 3 f ^^(x;P B My,P B )e- i ^ x °- Y ^e- iP ^-^x 
6(min[(xi)o, (x 2 )o\ - max[(2/i) , (2/2)0]) • 



here Pb = (wb, P) and ub = sJP 2 + m B ; m^-being the bound state mass 
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The argument of G function can be transformed as 

min[(xi) , (2:2)0] - max[(j/i) , (2/2)0] = 

min[X + rj 2 x , X - r/i£ ] - max[Y + r/22/o, *o - »7i2/o] = 

' X - lo - r/i£o - Wo , if xq > 0, 2/0 > ' 

X -Y - r)iXo + 7712/0 , if x > 0, 2/0 < 

X - Y + 7/2X0 - 7/22/0 , if x < 0, 2/0 > 

>X -Y + 7/2X0 + 7/12/0 , if x < 0, 2/0 < , 

X -Yq- ^\x \ - ^\y \ + ^(7/2 - 7/1) (x - 2/0) 
Besides, let us use its integral representation 



00 

e ( z) = ± J dk 



-ikz 



—00 



k + ie 



in which the variable change k — > Po ~~ is made. Then 



e(x - y - hx \ - hyo\ + ^(7/2 - m)(xo - 2/0)) = ^z j -^— Pn 



2tt P q -u b + ie 



x 



exp{-i(P - wb)[X - Y - ^\x \ - i|y | + 7^2 - 771) (ar - 2/o)]} 



inserting this into B, we get 



B(x, y-X-Y)= 7(2tt)- 4 / dP $(x; P B )$(y; P B ) 



g-iPpr-r) 



2w s (P - cjg + ze) 



x 



1. 1, 1 



exp{-7(P - uj b )[--\x \ - -|2/o| + -(7/2 - m)(xo ~ 2/o)]} = 



*(p; P B ) *(g; Pb) 

u B + ie) 



x 



J 2uj b \Pq — 

exp{-i(P -u B )[-^\x \ -^|2/o| + ^(r/2-m)(^o-2/o)]} , 

where in the reduced B-S amplitudes the transition to the momentum space 
was done: 



and 



$(x,P B ) = (2tt)- 4 / dpe- ipx $(p,P B ) 



<%, P B ) = (2tt)- 4 / dq e^%, P B ) 
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If we designate for brevity 

<S>(p-P B )<h(q-P B ) 



A( Pj q; P) = 



2cu B (P -u b + ie) ' 



f(x , y ) = ^(772 - r/i)(x - y ) - ^\x \ - ^\y \ 
then the momentum transform of B will be 

B(p,q,P) = (27r)- 4 | e ipx e~ iqy e iPX B(x, y; X)dx dydX = 

i(2tt)- 16 J dx dy dX dP' dp' dq'e l ^ p ' )x e- l{q - q ' )y x 

e i(p-nx A{p ^ q >. n exp { _ l{P , _ j B)f{x ^ yo)} = 

i(2tt)- 12 J dx dy dp' dq'e lip - p ' )x e- t{q - q ' )y A(p', q'; P') x 

exp {-i{P - u B )f{x , y )} ~ z(27r)~ 4 A(p, q; P) when P -> u B . 
Therefore, near the point Pq = u B 

j (27r)" 2 WB (P - WB + it) ' 

that is two-particle bound state contribution in G(p, q; P) has a pole at Pq = 
u B . Other |n > intermediate states will give a regular contribution at this 
point, if their masses differ from m B (and we assume that this is the case). 
So 

G{p,q;P) — ► — - 1 - — — — , whenPo -> w B . 12 

(27r) 4 2cj jB (P + ze) 

Bethe-Salpeter equation 

Let us take p ^ q, P — > in the equation (9) for Green's function. Taking 
into account (12), we get 

(27t) a 2lj b (P -u b + it) 

f A ,j t / p N * ^(^P^^^Pj?) 

•> (27T) 4 2u b (P -uj b + 16) 

This implies the Bethe-Salpeter equation for the B-S amplitude 

[MviPb + p) A 2 ( m P B - pT'Hp; P b ) = J dq' /(p, </; P B )<%'; P B ) . (13) 
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Normalization condition 



The normalization of the B-S amplitude can not be determined from the 
homogeneous B-S equation. It can be obtained in such a way. From (10) 
G • (K — I) = 1. Differentiating with respect to Pq, we find 



dG 

m 



\k-i) + g( 



dK dl 



that is 



dG 



= -G 



dP dP Q 
dK dl \ 



= 0, 



G 



dP Q \dP dP Q ; 

Let us consider this equation near the point Pq = ujb, where Eq.12 for Green's 
function is valid: 



i $(p; P B )Hq; Pb) 



{27ry2u B {P -u B + ie) 2 

d 



= — j dq dq 



dP 



(K-I) 



„ i <S>( P] P B Mq';P B ) 
{2ny2u B {P {) -u B + ie) 

i <S>{q"-P B )<h{q-P B ) 



x 



(2ny2u B (P Q -u B + ie) 



which gives the following normalization condition 



(2*)' 
or symbolically 



dp dq$(p; P B ) 



_8_ 

dP ( 



(K-I) 



J p =u B 



{p,qMq;P B )=2u; B = 2(P B ) 



i - (dK dl 

7$ 



(2tt) 4 \dP dPoJp^ 



$ = 2u B 



(14) 



B-S equation in the ladder approximation 

Let us consider the B-S equation in the ladder approximation for two scalar 
particles interacting via scalar quantum exchange. The interaction kernel 
in this case can be calculated according to standard Feynman rules: every 
vertex contributes ig and every scalar propagator 



A(X) = (2*)-* J - f P \ dp 

J —p l — l€\ 



Thus [A 



I(x, y; X — Y) = -gig 2 A (a:) S(x 1 - x 2 ) 5(x 2 - y 2 ) = 
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-g x g 2 AO) 5{X — Y + m (x - y)) 5(X — Y — 771(2; - y)) . 
Let us note that 

8{X-Y + m (x - y)) 5(X-Y- m (x - y)) = 

5[rn(x -y) + rj 2 (x - y)] S[X — Y — m (x - y)] = 5(x - y) 5{X - Y) , 
that is 

I(x, y-X-Y) = - gi g 2 A(x) 6(x - y) 8(X - Y) , 
its momentum space image being 

I(p,q;P) = (2tt)- 4 / e ipx e~ iqy e~ iPX I {x , y; X) dx dy dX = 
-g l92 {2n)- A J dx e^-^A(x), 

or 

I( P , q ;P) = - M -^ (p - q) = gg /l2 _ (p ^ )2 _ K • d5) 

Let us substitute this into (13) and also let us take Pb = 0. As a result, we 
get the rest frame B-S equation in the ladder approximation 

[mi+p 2 - {mPo + Po) 2 }[m 2 2 + p 2 - {r] 2 Po - p ) 2 ]<5>{p; P ) = 

A f <$>{q; P ) gig 2 

d( l-o 71 ^o—i A =T^- ( 16 ) 



itt 2 J fi 2 — (p — q) 2 — ie ' 167T 2 

This integral equation has a singular kernel and so standard mathematical 
tools are inapplicable to it. But the singularity in 

1 



v 2 - (po - qo) 2 + (p- q) 2 - ^ 

disappears if we suppose that p§ and qo are pure imaginary. The procedure of 
transition to such po, qo is called "Wick rotation" and it assumes two things: 
the analytic continuation of Pb) in the complex j^o-plane and a rotation of 
the dqo integration contour from real to imaginary axis. During this rotation 
the integration contour, of course, should not hook on singularities of the 
integrand. So it is needed to study the analytical properties of the Pb) 
amplitude. 
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Analytical properties of the B-S amplitude 

Let us consider the reduced B-S amplitude 

(2tt)- 3/2 $(x;P) =< O|T0i(?72z)02(-^)|£ >= 



where 
and 



e(x )f(x;P)+e(-x )g(x;P), 

f(x;P) =< Q\T<f>!(ri2p)<h(-mx)\B > 



g(x; P) =< QlTfai-mx^MlB > . 
Each of them can be transformed by using a complete set of states 

1 = H \ n >< n \ = E / dp\p, a >< p, a\ 



(a represents discrete quantum numbers of the state \n > and p is its 4- 
momentum). Taking into account identities 

fafox) = e 1 ^ >i(0)e-^ and <fo(-rnx) = e~ imxP (f) 2 (0)e imxP , 

we get 

f(x; P) = E I dp< 0\(f)i(r) 2 x)\p, a >< \p, a\4> 2 {-T]ix)\B >= 

a J 

E / dpe-^- 11 ^ < O|0i(O)b,a >< \p, a\(f) 2 {0)\B > , 
g(x;P) = J2[dp< Q\(fo(-riix)\p,ot >< \p, oi\<l>i(rbx)\B >= 
E / dpe- t{mP - p)x < 0\<p 2 (0)\p,a >< \p, a|<£i(0)|£ > , 



or 



f(x;P) = J dpe'^-^fip-P) , g(x;P) = J dpe- l{mP - p)x g(p; P), 
where we have designated 

f(p;P) = E < 0\M0)\P,<*>< \p,a\<h(0)\B > , 

a 

9(p;P) = E < O|0 2 (O)|p,a >< \p, a\M0)\B > . 



a 



The matrix element < O|0i(O)|p, a > is different from zero only then 
\p, a > has the same quantum numbers as the <p\ field quantum (otherwise 
0i(O)|p, a > will not have the vacuum quantum numbers). But among states 
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with quantum numbers of the first particle just this particle should have the 
smallest invariant mass, unless it will be not stable. So f(p; P) ^ only then 
Po > \Jm\ + p 2 . Therefore 

fix; P) = / dpe-^-^ x f(p; P)G(p Q - ^Jm 2 + p 2 ) = 

J dqe- iqx f( m P + q; P)G (q + m Po ~ V ™? + (Q + mP) 2 ) ■ 
This last equation, if we designate 

f(q; P) = ffaP + q; P) , u+ = \jm\ + (q + m P) 2 - m P , 
can be rewritten as 

fix; P)= j dq£° + dq e-^f(q; P). 

Let us note 

mi 

w+ > mi - 7/1 P = ; (rai + ^2 - Po) 

mi + m2 

In the rest frame Po equals to the bound state mass, therefore mi + m2 > Po- 
That is u + > and only positive frequencies contribute in fix; P). 

Analogously, < O|02(O)|p, a >^ only then \p,a > has quantum num- 
bers of the second particle and stability of this particle means g(p; P) = 
gip; P)O(p - sjml+p 2 ). So 

gix; P)= J dpe- l ^ 2p - p)x gip; P)O(p - V™2 + p 2 ) = 

J dqe- iqx gir] 2 P - q; P)G (r7 2 P - qo - ^Jrn 2 2 + foP - q) 2 ^ . 

If we designate giq; P) = g(7] 2 P - q;P), U- = t] 2 Pq - \jml + feP ~ Q) 2 ^ 
then 

gix;P) = / dqj^dq e- iqx giq; P). 

U- < r]2Po — m2 = m ^ m2 (Po — mi — m,2 ) < 0, so only negative frequencies 
contribute in gix; P) (in the bound system rest frame). 
As a result we get 

$(x;P) = (27r) 3 / 2 e(x ) / dqT dqoe-^fiq; P) + 
(27r) 3 / 2 6(-xo) / dqT {m dq e-^ x g(q; P). (17) 

J J — oo 

In the P = rest frame u + > 0, cj_ < 0, Therefore Eq.17 shows that from 
the xq > positive half-line $(x; P) can be analytically continued in the 
bottom half-plane (just then we will have falling exponent), and from the 
xq < negative half-line - in the upper half-plane. 

13 



Analytical properties of the B-S amplitude 
in the momentum space 

To deal with Eq.16, we need analytical properties of the Po) ampli- 
tude. Let us consider therefore the momentum space amplitude &(p;P) = 
/ dxe ipx §(x] P). Because of Eq.17 we have 

*(p; P) = 

(2tt) 3/2 J dqj™ dq J dx e'^-^ J™ dx e l{p °- q " )x " f(q; P) + 



(2tt) 3/2 / dq [ dq [ dxe~ i{ P-^ x [ dxoe l{p ^ qo)x °g(q; P). 

J J—oo J J—oo 

Integrals over dx, dq and dxo can be calculated, which gives 



z(2tT 



oo 



$(p;P) = 

f(Qo,p;P) 



dqo 

po - qo + ie 



(18) 



Po ~ Qo ~ ^ 

While integrating over dxo the following definition of the integrals was applied 

| o °° dx e i{ P°- q ^ Xo = | o °° dxoe l{p °- qo+ie)x ° , 

f° dxoe l{p °- qo)x ° ee f° dxoe l{p °- q °- ie)x ° . 

J—oo J—oo 

If now we try to continue &(p; P) analytically in the complex j^o-plane, every- 
thing will be O.K. except po — qo ± case. To avoid these points, cuts should 
be assumed in the j^o-plane. In the rest frame we will have the following 
picture: 



* 



Figure 5: Cuts in the p -p\&ne. 



In the remaining j9o _ P iane *&{p'i P) will be an analytical function. 
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Wick rotation 



oo 



-00 



In the r.h.s. of Eq.16 we have the following integral over dqo 

dqQ Hq;Po) 

v 2 - {po - qo) 2 + (p- q) 2 - ^ ' 

Its integrand has singularities shown in Fig. 6. 



it 



% 



Po > 

Figure 6: Integrand singularities 



Therefore 



/ dq c 



c - {po - qo) 2 + (p - </) 2 - «e 

where the integration contour is shown in Fig. 7. 



= 0, 





m > 



Figure 7: Integration contour. 
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The contributions from the two infinite quarter circles tends to zero. So 
Eq.19 means that the dqo-mtegidl over the real axis can be replaced by that 
over the C\ contour, which follows the imaginary axis and goes around a 
Po — \J/i 2 + (p — q) 2 + ie pole (if po > 0), or po + fi 2 + (p — q) 2 — ie pole (if 
po < 0). 
Therefore 



[m 2 +p 2 - (i^Po + p ) 2 ][m 2 + p 2 - ( m P -po) 2 }<S>(p;P ) = 
/ dq / aqo^ -; 

ITT* J J n* — Inn — 

Ci 



y 2 - {po - qo) 2 + (p- q) 2 - ^ 



(20) 



If now we begin to rotate po counterclockwise, both sides of Eq.20 remain 
well defined. When we reach the imaginary axis po = ip4, the disposition of 
the integrand singularities with regard to the C\ contour will be such: 




1 



Si 







It 



j!i 



ft > ® 



Figure 8: Disposition of the integrand singularities. 



As we see, the dangerous pole has left the C\ contour loop, so this contour 
can be straightened and it will coincide the imaginary axis. Therefore Eq.20 
becomes: 

[m\ + p 2 - {r)iP + ip 4 ) 2 } [ml + p 2 - {r] 2 Po - ipi) 2 ]^{ip4,p; P ) = 



J dq J dq - 



m 2 J _/ oo fi 2 - (ip A - q ) 2 + (p- q) 2 - ie 
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Let us make qo — > iq^ variable change in the integral over dqo: 

[ml + p 2 + O4 - imPo) 2 } [ml + p 2 + {pa + imPo) 2 ]^{m,p\ Po) = 
a f^7, Hm,q;Po) 



— I dql dq A ^gjjjj P o) 

7T 2 ' ^ V 2 + (P4-qA) 2 + (p-q) 2 ' 



If now introduce Euclidean 4- vectors p = (^,^4), g = ((f, ^4) and an amplitude 
Po) = ^{ipAiP'i Po) (which is defined by the analytical continuation of the 
B-S amplitude), then the above equation can be rewritten as 

[ml +p 2 + (pa ~ imPo) 2 } [m\ + p 2 + (pa + imPofMP; Po) = 
A mPo) 



7T 2 / 



(21) 



H 2 + (p - </) 2 ' 

This is the Wick rotated, rest frame, ladder Bethe-Salpeter equation. 
Salpeter equation 

One of the main peculiarities, which distinguishes the B-S equation from the 
Schrodinger equation, is the presence of the relative time xq (in the momen- 
tum space the corresponding quantity is the relative energy po). It is clear 
that one of the physical sources for its appearance is the retardation of the 
interaction. If there are no other physical grounds behind the relative time, 
its effect should disappear in the instantaneous approximation, when the re- 
tardation of the interaction is neglected. Let us check this for an example of 
two interacting spinorial particles, for which the B-S equation looks like 

-mi- 7(l)(»7i^ H 2 ~ rf^P^ ~ Pud] ^fo p ) = 

JdqI(p,q;PMq;P). (22) 

in the l.h.s. instead of full propagators we have taken free ones, that is 
A(p) = -^y. Besides #(p;P) = I dxe ipx ^(x; P) : and V(x;P) is a 16- 
component reduced B-S amplitude 

* a p(x; P) = (2irf 2 < 0\Tif> la (<n2x)il>2fi(-mx)\B > . 

7(1) matrices act on the first spinor index and 7( 2 ) matrices on the second 
one. 

The instantaneous approximation means that I(p, q; P) = I(p, q; P) does 
not depend on relative energy. Indeed, if this is the case, then 

I(x,y;X) = (2tt)- 8 | dp dq dPe- ipx e iqy e- lPX I(p,q-P) = 
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{27T)- 6 6{x )d{y ) J dp dq dPe l ^e- l ^e- lPX I(p,q-P) . 

Let us multiply both sides of Eq.22 over 7(°i)7(°2) and designate / = — 7(\)7(2)-f > 
then the equation in the rest frame P^ = (E, 0) becomes: 

[#i(p) - 771^ -Po][^ 2 (-p) - V2E + po]*(p; P) = 
JdqI(p,q;PMq;P), (23) 

where H(p) = a • p + /3m is a conventional Dirac Hamiltonian. The right 
hand side of this equation can be rewritten as 

/ dq~I(p, q- P)*(q- P) = J dqf(p, q; P)$(g; P) , 



where 



P) = J dq ^{q; P) = J dq dxe iqx ^(x; P) = 



-00 



2tt J dx5(xo)e-^(x , x; P) = 2tt j dxe-^(0, x; P). 

So $((f; P) is determined by a simultaneous B-S amplitude. 

To express the l.h.s. of the Eq.23 also in terms of $(<f; P), the following 
trick can be used. The projection operators 



A ± (p) = - 1± 



\A^ 2 + p 2 



have a property A±(p)H(p) = ±\/m 2 + p 2 A ± (p). using this, Eq.23 can be 
replaced by a system 

\jm\ + p 2 - T}\E — po — ie \jm\ + p 2 - ij 2 E + p - ie ^ ++ (p; P) = 
A { +\p)A?\-p) J dqI(p,q;P)mP) , 

^/m 2 + p 2 - 7/1 £" - po - ie -\jm\ + p 2 - ij 2 E + p + ^ + _(p;P) = 
A^A^-p) / dqf(p,q;P)mP) , 

-\/m 2 + p 2 - 7/1 £" - p + ie \fm%+ffi - r] 2 E + p - ie \P_+(p; P) = 
A^(p)Ai 2) (-p) / dqI{p,q;P)<5>{q;P) , 

~VW + p 2 - 7/1 £" - po + ie + P 2 ~~ ^P + po + ie * (p; P) = 

aL 1} (p)aL 2) (-p) y d#f(p, g; P)$(g; P) . (23') 
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Here ^ ++ (p;P) = A+\p)A+\— p)ty(p; P) and so on. Note the substitution 
m — > m — ie in the Feynman propagators (positive frequencies propagate 
forward in time, while negative frequencies - backward). 
By means of residue theory we get 

oo , 



m 2 + P 2 ~ V2E + po — ie 



-00 



-2m 



dpo \Jm\ + p 2 — r)\E — po — ie 

E — \fm 2j rp 2 — \fm 2j rp 2 
J dpo —\lm\ + p 2 - r\\E - p + ie —y/m% +p 2 - rj 2 E + p + ie 

2m 



-1 



00 



—00 



-1 



E + ym 2 + p 2 + \Jrri2 + p 2 



y dj9 \fm 2 + p 2 - r\\E - p - ie -\lm\+p 2 - T] 2 E + po + ie 



—00 
00 



/ g?77 +p 2 - r\\E - p + ie 



'm 2 , + p 2 — r]2E + po — ie 
Therefore for the $(<f; P) amplitude the following system holds 



= 0. 



E — ym 2j r p 2 — ym 2j r p 2 &++(p] P) = 
A^\p)A {2 \-p)[E - H^p) - H 2 (-p)Mp; P) = 
A { l\p)A { l\-p)— f dql(p,$P)mP), 

E + \jm\ +p 2 + yfml +p 2 $ (p; P) = 

A { l\p)A {2 \-p)[E - H^p) - H 2 (-p)]$(p; P) = 



-A { }\p)A^(-p)— I dql(p, q; P)m P), 

1 J 

$ + _(p;P)=$_ + (p;P)=0. 

To rewrite the last two equations in the same form, as the first ones, note 
that 



(2) 



2tt 



E — \Jm 2 + p 2 + ym2+ p 2 > , E + ym 2 + p 2 — ym 2 , + p 2 > . 
Indeed, if mi = 1112, these inequalities are obvious. Let m\ > 777,2, then only 

E — \Jm 2 +p 2 + \Jrri2+p 2 > 
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inequality needs a proof. But 




= m\ 



m 2 , 



\jm\+p 2 + ^Jml+p 2 mi + m 2 



and E > mi — m 2 is a stability condition for the first particle. Otherwise the 
first particle decay into the second antiparticle and \B > bound state will be 
energetically permitted. 

Because of the above mentioned inequalities, § + _(p] P) = §_ + (p] P) = 
equations are equivalent to 



K { l\p)K { l\-v) + A { +\p)A { -\-p) + A ( }\p)A { +\-p) + aL 1} (p)aL 2) (-p) = 1, 



So the relative energy is indeed excluded from the equation, but with the price 

of A ++ — A operator introduction. To understand why relative energy has 

left such a trace, it is useful to compare with the nonrelativistic case. 

Bethe-Salpeter equation in the nonrelativistic theory 

To derive the nonrelativistic Bethe-Salpeter equation, one can use the fact 
that quantum field theory in many respects is similar to a second quantized 
many particle theory. In graphical representation, to the second quantized 
Hamiltonian 




A^(p)A {2 \-p)[E - Hi(p) - H 2 {-p)\m P) , 



= E + \jm\ +p 2 - \jm\ +p 2 P) = 

A { l\p)A { l\-p)[E - Hi(p) - H 2 (-p)]<5>(p; P) . 
If we sum all four equations for P) and use 



then we get the Salpeter equation 






dx dy t)^ + (y,t)V(\x - y\)^(x, t)^(y,t) 
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there corresponds the free propagator 



A(x - y) = 



dp 



■ip(x—y) 



(2tt) 4 J p - (p 2 /2m) + ie 



and the pair interaction (instantaneous) with the potential V: 

-id{x - yo)V(\x-y\) . 

Everywhere in the derivation of the B-S equation , (j)(x) field operator can be 
replaced with the second quantized operator fy(x). As a result, we end with 
an equation for the following nonrelativistic Bethe-Salpeter amplitude 

$0; P) = (2tt) 3/2 < 0\T^/{r] 2 x)^(-r] 1 x)\B > . 

Namely, in the ladder approximation, the interaction operator is 

I(x h x 2 ; 2/1, y 2 ) = -iS(x )V(\x \)5{x 1 - yi)S(x 2 - 2/2) = 

-i6(x )V(\x\)6(x - y)6(X -Y) , 
with the corresponding expression in the momentum space 

/(p, q- P) = (2tt)- 4 J dx dy dXe ipx e~ iqy e iPX I(x, y; X) = 



-z(2tt)" 



-i(p-q)-x 



— I 



■ x V(\x\)dx=-^y 4 V(p-q) 



Besides, the momentum space free propagator looks like A(p) = i/(po — 
p 2 /2m J r ie), therefore the B-S equation, in the rest frame and ladder approx- 
imation, will take the form 



p 



(t]iE + po) - h ie 

2m 



(mE-po)-^ + ie 



$(p; E) = 



i(2n)- 4 fdqV(p-q)<P( q] E) 

(note that mi = m 2 — m and so 771 = 772 = - 

Let us introduce integrated over dpo amplitude $(<f; E) = J dqo $(q; E). 
For it the Salpeter equation holds 



-00 



m e) = 



dpc 



[rjiE + Po- p 2 / 2m + ie] [r) 2 E — po - p 2 / 2m + ie] 



x 
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(2*)' 



/ dqV(p-q)$(q;E) 



The square bracket integral can be evaluated via residue theory and it equals 
—2iri/(E — p 2 /2m). So we obtain the following equation 



E - 



p 



2m 



®{p;E) = (27r)- 3 | dqV{p-q)${q;E) 



(25) 



But this is just momentum space Schrodinger equation! Indeed, in config- 
uration space 



= / e -^' f |)(f) dx, V(p) = J e-&*V(\x\) dx , 



we will have 



-- + V(\x\) 
m 



!>(£) = E$(x) 



Thus, in the nonrelativistic theory, the relative time can be excluded without 
any trace. So its introduction is purely formal. 



Physical meaning of the relative time 

What is the crucial peculiarity, which distinguishes the above given nonrela- 
tivistic model from the instantaneous approximation of the relativistic one? 
It is the propagator! There are no antiparticles in the nonrelativistic case. 
Therefore the propagator describes only forward propagation in time, that 
is we have retarded Green's function: A(x) = 0, if xq < 0. This boundary 
condition demands po — > po + ie prescription for the propagator poles. Let 
us see what will be changed in the Salpeter equation derivation if we replace 
the Feynman propagator i(p + m)/(p 2 — m 2 + ie) by the retarded Green's 
function i(p + m)/[(po + ie) 2 — p 2 — m 2 }. Instead of Eq.23' system, we will 
have 

V ++ (p;P) = 



*+-(p; P) = 



*-+(p; P) = 
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VW + p 2 - rjiE - p - ie ^ml + p 2 - rj 2 E + po - ie 

A ( +\p)A ( +\-p) J dqI(p,q;P)mP) , 
\Jm 2 + p 2 - r\\E - po-ie -\jm\ + p 2 - r\ 2 E + p - it 

A^\p)A?\-p) J dqf(p,q;P)mP) , 
—\jm\ +p 2 - r\\E - po- 



ie 



m 2 +p 2 — r]2E + po — it 



A { }\p)A^\-p) J dql(p, q; P)$(q; P) 



ttt, 2 + p 2 — r\\ E — pq — ie — \jm\ + p 2 — r\iE + p§ — ie \I/ (p;P) = 



A { l\p)A {2 \-p) J dql(p, q; P)$(q; P) 



(the propagators, which appear in the B-S equation, have T]iP+p and 772P— 
as their 4-momenta. Therefore the transition to the retarded Green's function 
means the following replacements 

(771 P + p ) -> (771P0 + Po) + , (mPo - Po) -> (m p o - Po) + ^ 

that is ie has the same sign, as E). Integration over dpo can be performed 
using 

dp 



00 



/ 

J (a — Da — 



= 2ni 



-00 



1 

a + b 



(a-p - ie)(b + p - ie) 
and we get the system 

' $ ++ (p; P) ee A ++ [£ - ^(p) - H 2 (-p)]<S>(p; P) 
2n r 

dqI(p,q;P)<S>(q;P), 



'm\ + p 2 — 


^777-2 + p 2 




= A ++ 


'm\ + p 2 + 


y^m 2 + p 2 




= A + _ 


'm 2 + 7J 2 — 


\Jm2+p 2 




= A_+ 



2tt 

z 

z 



j dqI(p,q-P)$(q-P), 



JdqI(p,q;P)mP), 



E + v/m 2 + 7j 2 + Vmf+^j $_ p) = A __[p _ _ # 2 (_p)]$(£ p) 

z 



= - j dqUiiq-.PYHq-.P). 



Summing these four equations, we get the Breit equation 

2-7T r 

[E -HM- H 2 (-p)Mp; E) = — dql(p, q- P)$(q; E) 



(26) 



which is a direct generalization of the two particle Schrodinger equation 
Eq.25. (The nonrelativistic Hamiltonian H(p) = p 2 /2m is replaced by the 
Dirac Hamiltonian H(p) = a ■ p + [3m and —i(27r) 4 I(p,q) plays the role of 
potential). Once again, the relative time disappears without any trace left. 
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Therefore, the second (and more important) source for the essential rela- 
tive time dependence of the B-S amplitude is the existence of antiparticles, 
that is the possibility for particles to turn back in time and propagate back- 
ward. Because of retardation of the interaction, relative times of the order of 
bound system size will be significant, while the forward-backward motion in 
time makes essential configurations, for which individual routes in time are 
very different for bound state constituting particles, and so the relative time 

is large. A ++ — A operator in the Salpeter equation just corresponds to 

the contribution of these configurations to the bound state amplitude. 



Wick-Cutkosky model 

This model corresponds to the B-S equation in the ladder approximation for 
two scalar particles interacting via massless quanta exchange. After the Wick 
rotation, the corresponding rest frame B-S equation takes the form 

[ml +P 2 + (p 4 " iriiE) 2 } [m\ + p 2 + (p 4 + irj 2 E) 2 }<S>(p) = A J dq ^ (27) 

where E is the bound state mass and p, q are Euclidean 4- vectors. 

To investigate the mathematical structure of Eq.27, let us consider at first 
the simplest case mi = = m and E = (although, in the rest frame, 
E = is, of course, unphysical: massless bound state doesn't have rest 
frame). Eq.27 then becomes 

iP i +mrm= ±j dq _m_. (28) 

Let us show, that one of its solutions is (j)(p) = (p 2 + ra 2 )™ 3 . We have 

^-^(44)- 

dx r dx 

[B + {A-B)x] 2 = J [xA+{l-x)B] 2 ' 

A-n-D-m = [ L ) ° ° / A-l d-1\ 

^ (n - l)\(m - 1)! dA n - 1 OB™- 1 { ) ' 

and we get the following Feynman parameterization 

A-B-™ = ^ n+m - i y- } g^q-sr- 1 dx m 

{n-l)\(m-l)\{[xA + (l-x)B] n +™ K } 



besides 
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In particular 

3 
So 



[(p - qf}-V + rnT' 6 = 



2-1-3 



1 — x) 2 dx 



3(1 — x) 2 dx 



o 



[x(p — q) 2 + (1 — x)(q 2 + m 2 )] 4 ^ [(g — xj^) 2 + (1 — x)(m 2 + xj9 2 )] 4 ' 

dg g(g) r rf ? 3(1 -xfdx 

J I i(i xi)) 2 ■ (1 x)(m' 



If now we use 

A- 4 = 

then we get 



(p - q) 2 

1 d 3 /I s 
~3! 

l 



[(q — xp) 2 + (1 — x)(m 2 + xp 2 )] A 

i ,q3 00 -i oo 

u j e - aA da = -ja 3 e- aA da, 



3!6U 3 



3(1 - x) 2 dx 



[(q — xp) 2 + (1 — x)(m 2 + xp 2 )} 4 

1 oo 

- J {1 — x) 2 dx J a 3 da J dqexp {— a[(q — xp) 2 + (1 — x)(m 2 + xp 2 )}}. 



The Gaussian integral over dq equals 

J dqexp {— a[(q — xp) 2 + (1 — x)(m 2 + xp 2 )]} = 



7T 



a 1 



exp {— a(l — x)(m 2 + xp 2 )}, 



and we are left with an integral over da of the type 



OO o oo 

/ ae^da = -— / e^da = ,4~ 2 
OA J 



As a final result, we get 



dq- 



7T 



2 1 



7T 



(p — q) 2 2 ■> [m 2 + xp 2 } 2 2m 2 



[m +p) 



2\-l 



Inserting this into Eq.28, we will see that $>(p) = (p 2 + m 2 )~ 3 is indeed a 
solution and the corresponding eigenvalue is A = 2m 2 . 

The most interesting thing about this solution is that we can indicate its 
analog in the nonrelativistic hydrogen atom problem. 

The Schrodinger equation 

A 



2m 



+ V{r) 



#(f) = Em{r) 
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in the momentum space becomes an integral equation 



(p 2 - 2mE)^(p) = 



2m 



(2tt) 3 / 2 



dq V(p — q)^(q) , 



where V(q) = (2n)-^ 2 J e- l ^{r)dr and V(p) = (2n)- 3 / 2 1 e- l ^V{r)dr. 
A(r _1 ) = —An5{r) identity indicates that V{r) = —e 2 /r Coulomb potential 
has V{p) = — (47r/(27r) 3 / 2 )(e 2 /p 2 ) as its momentum space image. Therefore 
a nonrelativistic hydrogen atom is described by the equation: 



(30) 



where p 2 } = —2mE (we consider a discrete spectrum and therefore E < 0). 

One of the solutions of this equation is ^(p) = (p 2 + Pq) ~ 2 • Indeed, in the 
similar way as above we get, after integrating over dq 

Jdq[(p-q) 2 ]- 1 [p 2 +p 2 }- 2 = 

1 oo 

^3/2 J ^ — x) dx J v^exp {—a(l — x){p\ + xp 2 )} da . 
o o 
Integral over da can be evaluated in such a way 



/ Va~e~ aA da = 2 J t 2 e~ eA dt = -2^- J e 



OA 



- t2A dt = 



d 



'7T 



1 [W 



therefore 



f Jll 
J (r>- 



7T 



dA\ A 2A\A' 

dx 



2 1 



(p- q)' 



-2tt< 



d{pl) 2 H i-x dx ~ n d( P 2 yJ\ i-x p dx 



2 I {pi + xp 2 )^{l - x){pl + xp 2 ) 

~ d 2 _ 

ha 





1*2 A rf 



pfL oo , , 



7T 

~2 



2 oo 



7T 



2(p 2 Q +p 2 ) 



oo 

J x~ 3/2 dx = 



pl + P^t-p 2 ]-*' 2 dt = 



7T 



Pi 



Po(Po + P 2 ) 
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Inserting in Eq.30, we see that if me 2 /po = 1, then the equation is fulfilled. 
Thus ip(p) = (p 2 + Pq)~ 2 corresponds to the E = zz2 ?p- ground state of 
hydrogen atom. 

If the analogy between the found solutions of Eq.28 and Eq.30 is not ac- 
cidental, one can expect that the same methods, which are used in dealing 
with hydrogen atom, will be useful also for Eq.28 and maybe even for Eq.27. 
In particular, it is well known that the nonrelativistic hydrogen atom pos- 
sesses a hidden symmetry and its study is more easily performed in the Fock 
space, where this symmetry becomes explicit. For the sake of simplicity, let 
us illustrate the Fock's method on an example of a 2-dimensional hydrogen 
atom. 



Fock's method for a 2-dimensional hydrogen atom 

The Scrodinger equation in the 2-dimensional momentum space looks like 

{p 2 - 2mE)^{p) = -— f dqV{p- q)^{q) . 
We have the following connection with the configuration space 

#(p) = _L f e - if? ^{r)dr, V(p) = — — / e-® r V(r)dr. 

For the Coulomb potential V{r) = —e 2 /r the above given integral diverges. 
But let us note that as bound state wave function is concentrated in a finite 
domain of space it should not feel a difference between Coulomb potential 
and a (—e 2 /r)e~ ar potential for sufficiently small a. Therefore, at least for 
bound states, the momentum space image of the Coulomb potential can be 
defined as 

V{p) = lim — J e-® ? [--) e- ar df = 



3 2 oo 2tt 

27T 



-—lim / e- ar dr I e^ rcose d6 



Using 



2tt 



/ e- iprcose dO = 27TJ {pr) , / e- ar J {pr)dr = 



we get 



e 2 



v ' y/p* + a 2 \p\ 
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So the Schrodinger equation for the 2-dimensional Coulomb problem, in the 
case of discrete spectrum E < 0, will be (pi = —2mE): 

(f + vlmP)~jP^dq. (31) 

7T J \p-q\ 

The 2-dimensional momentum space can be mapped onto the surface of the 
3-dimensional sphere (which we call the Fock space) by means of the stereo- 
graphic projection. The stereographic projection transforms a 2-dimensional 
momentum p = p x i + p y j into a point on the surface of the 3-dimensional 
sphere where this surface is crossed by the line which connects the south pole 
of the sphere with the {p x ,P y ) point in the equatorial plane. This is shown 
schematically in Fig. 9 below. 



4 




Figure 9: The stereographic projection. 
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Let the radius of the sphere be po- If the polar coordinates of the vector 
p are (p, 0) when the point P will have the spherical coordinates (po, G, 0), 
where the angle G, as it is clear from the Fig. 9, is determined by equation 
|p| = Po ^5 f ( ( 9/2). Therefore 

G 2 p 2 — p 2 

cos G = 2 cos 2 — — 1 = ^7— 7-7 — 1 = » 

2 l + ^ 2 (e/2) pg + P 2 

and 

• n /1 2?^ 2 PoP 

sin G = v 1 — cos 2 G = 



Po +P 

As for the Cartesian coordinates P x = po sin G cos 0, P y = po sin G sin 0, P z = 
Po cos G of the point P, we will have 

2p 2 p cc 2p 2 p y _ pg-p 2 

£ 2 1 9 ' y 2 1 9 ' z 2 1 9-^0 • V oz 7 
Po + P Po + P Po+P 

Let P and Q be stereographic projections of the vectors p and <f, and let a be 
an angle between 3-dimensional vectors P and Q. When, as |P| = \Q\ = po, 
the distance between P and Q points in the 3-dimensional space will be 
\Jpo + Po — 2po cos a = 2po sin (a/2). On the other hand, the square of this 
distance equals (P x — Q x ) 2 + (P y — Q y ) 2 + (P z — Q z ) 2 , therefore 



2p oS in-l 2 = [ 2pZ ° Px 2pZ ° qx V I [ 2pZ ° Py 2p2 ° qy 1 -+ 

v 2/ Ipo + p 2 Po + q 2 ) \pI + p 2 Po + q 2 




2 



(Po + P 2 ) (Po + 9 2 ) l Po + P 2 Po + 9 s 

/ Po + ff 2 ^ pg+P 2 _ 2 pV 
Po Po + P 2 Po Po + q 2 Po 



4pp P 2 ( 2 , 2N , Q ( 2 , 2n - PV 

/ 2 , 2W 2 , 2\ 1 ~2 (PO + ^ ) + — PO + P ) - 2 P * q ~ 

(po + p 2 )(Po + r) [po Po p 2 

(p- 9)' 



4P0 1 



(Po + P 2 )(Po + 9 2 ) 



But 



4p4 1 2 e 2 e' 

= cos — COS 



(p§ + p 2 )(p§ + ^ 2 ) (i + V(e/2))(i + v(eV2)) 2 2 
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where 0' spherical coordinate corresponds to the point Q. So we get 



6 9' 



a 



\P~ q\=Po sec — sec — sin 



(33) 



Now let us express dp = dp x dp y through the solid angle element dQ, = 
sin 6 d@ 



dQ = — d(cosO) dcf) = —d 



p\~P T 



4^o 



(rf + F 2 ) 2 



p dp d(f) = 



\ 2po ' 


dp = 


[2 2 
— cos 


0- 


-Po + P 2 - 




2". 



n2 



Thus 



= — sec 
4 2 



4@ 



(34) 



With the help of Eq.33 and Eq.34, it is possible to rewrite the Schrodinger 
equation (31) in the Fock space 







me 2 f ^{q)sec 3 e'/2 



sec 3 ^(p) = ^f— I IMl^l dQ' 
2 yFJ Aitp Q J sin a/2 



or, if a new unknown function $(£1) = sec 3 %^>{p) is introduced 



*(n) = — / 



47T]9o ^ sin a/2 



(35) 



This is the Fock space Schrodinger equation. It is invariant with respect to 
a 3-dimensional rotational group 50(3) (the angle a measures an angular 
distance between two points on the sphere surface and so it is unchanged 
under rotations of the sphere). In the momentum space this symmetry is 
hidden: Eq. 31 is explicitly invariant only under a smaller group 50(2). 

A solution of Eq. 35 is proportional to the spherical function Y\ m . Indeed, 
from the generating function of the Legendre polynomials 



oo 



(1 - 2ecosa + e 2 )~ i/2 = £ e'^(cosa) 

1=0 



we will have for e = 1 
1 



oo 



47T 
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and therefore 



Yan) = 2 ^J^-dn'. (36) 

Sir J sin (a/2) 

So the general solution of Eq.35 has a form <&(£!) = J2 l m ^i a m Yi m (Q) , if the 
eigenvalue condition (me 2 /po) = / + (1/2) is fulfilled. The corresponding 
energy eigenvalues are Ei = -(fjfp. Each of the energy levels is (21 + l)-fold 
degenerate, and this degeneration is of course caused by the above mentioned 
hidden symmetry, which in the Fock space becomes explicit. 



Stereographic projection (general case) 

To generalize the Fock's method for equations (28) and (30), one needs first 
of all a definition of the stereographic projection when the momentum space 
dimension is greater than two. But the above given definition of the stereo- 
graphic projection admits in fact a trivial generalization for the n-dimensional 
momentum space. 

Let a momentum space point p = (pi, . . . ,p n ) is transformed into a P = 
(Pi, . . . , P n +\) point on the sphere surface by the stereographic projection. 
sp and pP vectors are parallel because P lies on the line connecting the 
point p to the "south pole" s = (0, . . . , 0, — po) of the sphere. Therefore 
(Pi — pi, . . . , P n — Pn-, Pn+i) = k(pi, . . . ,p n ,Po) ? where k is some constant. So 
Pi = (k + l)pi, i = 1 -T- n and P n+ i = kpo. But the point P is on the sphere 
surface and therefore Pf + • • • + P^ + P^ +1 = Pq. This enables us to determine 
k from the equation (k + l) 2 p 2 + k 2 pl = p 2 } as k — p °~ p 

Therefore a generalization of the Eq.32 is 



Po+P 2 ' 



_ 2p 2 Pl 2p 2 p n _ pj-p 2 

Po + p z Pq + p z Po + p z 

It is also straightforward to generalize Eq.33 

W- q\ = Po sec ^ sec ^ sin ^ . (38) 

(p and q are n-dimensional vectors, a is the angle between P and Q (n+ 1)- 
dimensional vectors). A proof of Eq.38 looks like this: 

2p2 s in 2 1 = p 2 (l - cosa) = p\ - P ■ Q = 
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Pi 



W' q 



{pi + p 2 ) (p§ + g 2 ) (p§ + p 2 ) {pi + g 2 ) 

e 



Oo - P 2 ) (Po - 9 2 ) 2 _ 2 rfb 2 + 9 2 - 2p • g) 
Po - 



Po 



Po 



Po + P 2 Po + 9 2 v " J/ 2 2 

Now it is necessary to generalize Eq.34. A connection between spherical and 
Cartesian coordinates in the n-dimensional space reads 



:{p — q) = 2 cos z 



{pl+p 2 ){pl + q 2 ) 

(p- q) : 



(-y 



n 


= r sin O n - 


-i sin6 n _ 


-2 • 


.sin 63 sin 62 cos 6i 


^2 


= r sin G n _ 


-isin6 n - 


-2 • 


.sin 63 sin 62 sin 61 




= r sin 6 n _ 


-isinG^ 


-2 • 


. sin 63 cos 62 , 


r 4 


= r sin n _ 


.isin6 n _ 


-2 • 


. sin 64 cos 63 , 




= r cos 6 n _ 


-l • 







At that for the n-dimensional volume element S n ^r '= dri . . . dr n we will have 



d(")f = r n ~ l dr d^Q = 



r n 1 sin n 2 6 n _i sin n 3 6 n _ 2 ... sin 6 2 dQ\ . . . o?6 n _i 



(39) 



This equation can be proved by induction. As for n = 3 it is correct, it is suffi- 
cient to show that from its correctness there follows an analogous relation for 
the (n + l)-dimensional case. Let us designate r = yV 2 + . . . + r 2 +1 and p = 



r\ + . . . + r 2 = rsin6 n , when d {jl+V) f = d^r dr n+ \ = p n l dp dr n+ \ d^Q, 
but 



cos G n sin G n 
— rsinO n rcos6 n 



dr dOn = r dr dO 



n 5 



therefore 

= (rsine n ) n -Vdrde n d (n) Q = r n dr sin 71 " 1 O n dO n d {n) n = r n drd {n+1) n 

At a stereographic projection 

2p p 



sinO„ = 



Po + P z 



and cosG„ = 



2 2 
Po +P 2 



therefore 

d (n+ ^n = sin"- 1 e n de n d (n) Q = - sin n ~ 2 e n d(cosO n )d (n) n = 



' ZpoP 

kpI + p 1 



n-2 



d 



Pp-P r 
vP0 + P 2 
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Po \Po + P J \PoJ 2 

So, Eq.34 is generalized in such a way 

dF>p= (^-Jsec 2n ^S n+ ^n . (40) 

By using Eq.38 and Eq.40 it is possible to rewrite equations (28) and (30) 
in the Fock representation. In particular, Eq.30 for a hydrogen atom, after 
an introduction of the new $(£1) = sec 4 ^.^(p) unknown function, takes the 
form 

me 2 r W) 



<p(n) = — — / , 2 \ J — dQ! (41) 
v 7 8^0^ sin 2 (a/2) v ' 

where dtt = d^Q = sin 2 6 3 sm02d0 1 d0 2 d0 3 . 

As for the Eq.28, it can be rewritten, for the new $(p) = sec 6 ^$(p) 
function, as 

= — ~o~~9 I x dn' (42) 

V ; 16ir 2 p 2 J sin 2 (a/2) V ; 

here po = m and dQ — d^Q = sin 3 64 sin 2 63 sin 02<i0io?02(i03(i04. Note 
that, as it turned out, Eq.28 possesses a hidden 50(5) symmetry, which 
became explicit in the Fock space. 

So, as we see,there is indeed a close analogy with the nonrelativistic hy- 
drogen atom problem. 

To solve equations (41) and (42), one needs a generalization of the Y\ m 
spherical functions for the n > 3 dimensions, as the example of Eq.35 sug- 
gests. 



Spherical functions in a general case 

It is well known that in a 3-dimensional space spherical functions are con- 
nected to solutions of the Laplace equation. Namely, yi m (r) = r l Yi m (Q,4>) 
harmonic polynomials (solid harmonics) are homogeneous polynomials of 
rank / with regard to the variables r x ,r y ,r z and obey the Laplace equation 
^yimif) = 0. For the multidimensional generalization we will use just this 
property of the spherical functions. 
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Let us define n-dimensional spherical function fi (0 n _i, . . . , Gi) by 

requirement that 3^™_\ i x (r) = r^Y^ fl (G n _i, . . . , Gi) turns up to be a 
harmonic polynomial, that is obeys the equation 



A ( " ) y£... Il (r-)=0, 



(43) 



where = + . . . + is n-dimensional Laplacian. For it we have the 
following decomposition into a radial and angular parts: 



A («) = £ + _ Lfin) 



(44) 



being the angular part of the Laplacian. 
Let us prove Eq.44 by induction. For n = 3 it is correct. So let Eq.44 is 
fulfilled and consider 



<9 2 



d 2 n — I d 1 



<9r 2 



n+l 



dp 1 



p dp p' 



d 2 



dr 2 



n+l 



here p = \jr\ + . . . + r 2 . From r n+ i = r cos G n and p = r sin G n we get 



5 d dr ^ 



dp drdp dO n dp 

d d dr d dO 

+ 



d dG n 1 5 

= sin G n — + - cos G n — — , 

or r oB n 



d 



d 



dr n+1 dr dr n+ i dG n dr n+ i 
and for the second derivatives 

d 2 ■ 2r, d " 1 2r> d 2 
= sin G n — r H — -cos G 



= cos O n j- - - sin O n — — , 
c*r r oO n 



dr 2 



dG 2 



2 ^ ^ ^ 
-^cosG n smG n — -■ 
r 0\j n 



2 

r 



d 2 
n drdO 



1 2 d 

-cosQ„ sinQ„ - _ — h-cos G n — 

n r or 



d 2 d 2 1 d 2 2 5 

_ _ - = cos 2 G n — -z + -2 sin 2 O n — — + cos G n sin G n — 
or n+ i or 1 r l OOi r l OO 



n 



2 d 2 

- cos O n sin G n — — — 
r OroOn 



1 



5 



+ - sin 2 G n -- 
r or 



Substitution of these expressions into A n+1 gives 

j2 



A (n + 1) = ^1 + _ I 

^^^2 ^ (^T* y*2 



sin 



- (n - l)c*0 G n 



5 2 



2 G 



dO n d0 2 n 
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So the following recurrent relation holds 

S (n+1) = _ ( n _ ^ ' ' 



«9 2 



sin 2 6 n 



"de n del 



(45) 



From Eq.43 and the decomposition (44) it follows that a n-dimensional spher- 
ical function obeys an eigenvalue equation 

^11^(^-1, ■ ■ ■ , 6i) = In-^k-! + n - 2)y£ ) 1> ... i/l (e n _i, • • • , 6i). (46) 
Eq.45 suggests that the spherical function has a factorized appearance 

v£:;!(e B , ...,e,) = /(e„)y£,... A (e^i, . . . , e,) , 

Insertion of this into Eq.46 gives the following equation for /: 



^- 1 + n -^l- ( n-i)^e a 



6> 2 



sin 



2 e 



n 



/(Bn) = 



/„(i„ + n-l)/(e„) 



Since 



5 



©nT^" = - COS 6 n 



d 



d cos 6 n 




and 



«9 2 



d d 



(-sin6 n ) 



d cos 6 r 



-cos6 n - — -sin6 n - 

<7fc) n <7 cos B n 



— cos 6 



d cos 6 n 
8 



d cos 6 n 

this equation can be rewritten as 



+ sin 2 6 r 



s 



d(cos 6 n ) 2 



(l-x z 



d 2 f 



df 



dx 2 nX dx 



Iniln + n - 1) - 



(n-lftt-l + U-2) 
1 — X 2 



f = o 



where x = cos6 n . Let us take f(x) = (1 — x 2 ) ln - 1 l 2 g(x), when the equation 
for g will be 



(1 - x 2 ) 



d 2 g 



dg 



2 L n + 2/ n _i]x— + (l n - l n -i)(l n + /„_i + n - l)g = 



It should be compared with equation defining the Gegenbauer polynomials: 

(1 - x 2 )^C^(x) - (2a + l)x^C$\x) + N(N + 2a)C^(x) = . 
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As we see, g{x) is proportional to Cf n _i + ^ n l ^ 2 \x). Therefore a multidimen- 
sional spherical function can be defined through a recurrent relation 

i£+ 1 >(e n ,...,e 1 ) = 
a u ^ S in'»- e n ctt + ; n - 1)/2) (co S e„) vitL^e^!, . . . , ej , (47) 

where Ai n j n l constant is determined by the following normalization condition 
/ \Y^j^\ 2 d^ n+1 ^Q, = 1, which, if we take into account 

S n+1 ^n = sin 71 - 1 e n de n d {n) n = - sm n ~ 2 o n d(cose n )d {n) n 

can be rewritten as 

A 2 l j \c? n -? +{n - 1),2 \x)} 2 (1 - x 2 )^ n+2l -^dx = 1 . 
-1 



But for the Gegenbauer polynomials we have 

I (1 X > [Cn (X " dX ~ JV!(JV + a)[r(a)P ' 

therefore: 

,2 = gn ~ /n-i)!(2f„ + n - lWVn-1 + (n - l)/2)] 2 

This expression together with Eq.47 completely determines multidimensional 
spherical function. 

Namely, 4-dimensional spherical function looks like 

W©3, e 2 , Bi) = 



n + i)( n-i)wr ( sin e 3) -e + / ) (co S e 3 )y ( „ l( e 2 ,e 1 ) , (49) 

\ 7r(n + / + lj! 



and a 5-dimensional one: 



we* e* 6*9,) = ;»^ff I ^x 

(sin e 4 )"cK /2) (co S e 4 )y„ im (e 3 , 92,90 , (so) 

It is clear from these equations that I < n < N . In general, as a lower index 
of the Gegenbauer polynomial coincides to the polynomial rank and so is 
not negative, from Eq.47 it follows the following condition on the quantum 
numbers l n > l n -\. 
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For spherical functions the following addition theorem holds: 

e' - e e y£ ) 1 ... ll (e n _ 1) ...,e 1 )y 1 it i ,(e;_ 1) ...,e' 1 ) = 

*n-2=U /2=U/i = — / 2 

21 #r(¥)c ((n-2)/ 2)(cosa)) (51) 

where a is the angle between n-dimensional unit vectors, which are defined 
by spherical coordinates (61, . . . , O n _i) and (0' l5 . . . , 0^_i). 

If n = 3, Eq.51 gives the addition theorem for the 3-dimensional spherical 
functions Y/ TO , since cf (cosa) = P/(cosa). So let us try to prove Eq.51 by 
induction, that is suppose Eq.51 is correct and consider (n + l)-dimensional 
case: 



E y£ 1 1) (e„,...,e 1 )KS, 1) *(e; i ,...,e' 1 )= £ ^L.- x 



S i„'-> e„ sin 1 - e;cfcl + , ( "- 1)/2) (co S eocfc^'^'Ox* 9'„) x 



2"- 3 (2/ n + n-l)r((n-2)/2) 



E CL.J0.-i, ■ ■ ■ , ©i)&(©U •••,©;) = 

k-2,---,h 



47r (n/2+l) 



X 



l„-i=0 r ('« + 'n-1 + » - 1) 

[r(?n-i + (n - i)/2)] 2 S in'-' e„<^r^- 1)/2 W e„) x 
sin'- e;,c«:r ( ;! 1 ( "- 1)/2) (cose'„)c«»- 2)/2) (co SW ) . 

u is the angle between n-dimensional unit vectors and it is connected with the 
angle a between (n + l)-dimensional unit vectors as cos a = cos0 n cos0' n + 
sin Q n sin Q' n cos u (if e is (n + l)-dimensional unit vector, determined by 
(0i, . . . , n ) angles and / n-dimensional unit vector, determined by 

— * — * 

(01, . . . , n _i) angles, when cosa = e • e ' , coso; = / • /', and e = 

— * 

cos Q n e n +i + sin Q n f). 

Using the addition theorem for the Gegenbauer polynomials 

(cos 6 cos 0' + sin 6 sin 0' cos 0) = 

sin- e'C^r^cos e / )Ci- 1 /2)( CO s 0) , 
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we get: 

E y£, 1) (e„,...,e 1 )y£; 1) *(e;,,...,e' 1 ) = 

on _ 3 2/ n + n- l [r((n - l)/2)] 2 r((n - 2)/2) ((n _ 1)/2) 

2 47r («/2+i) f(^2) z « t COSQ 0' 

but T(z)r(z + 1/2) = ( 7 r 1 / 2 /2 2z - 1 )r(2z) formula gives T((n - l)/2)r((n - 
2)/2) = (7r 1 / 2 /2 n ^ 3 )r(n - 2) and we end with the desired result: 



lu- 



ll 



2Z„ + n — , , , „ .. ^((„-i)/2) 



^ 7 ^r((n-l)/2)Cr l " / '(co S a). 

In a 4-dimensional space the addition theorem for the spherical functions 
looks like 

n I 

E E K lm (e 3 ,e 2 ,e 1 )^(e' 3 ,e' 2) e' 1 ) = 

1=0 m=-l 



n+1 f1l/ , n+lsin(n+l)a 

— — C£ J (cosa) = — ; , 

2-k 1 2tt z sin a 

Using this, we can prove, that a solution of Eq.41 is proportional to a spherical 

function Y n i m (Q). Indeed, from the generating function of the Gegenbauer 

polynomials (1 — 2ecosa + e 2 )~^ = E^Lo C^\cosa)e n , we get in the e = v = 1 

case 1/4 sin 2 (a/2) = E^Lq C^(cosa) Note that it is not rigorously correct 

to take e = 1 because in this case the series stops to be convergent. But 

™ , ^sin(n+l)a sin ((JV + l)/2)a sin (Na/2) 

E Ci. ; (cosa) = E : = nm : —, — 77^ ■ 

u =q n=0 sin a n->oo sin asm (a/ 2) 

of course a limit doesn't exist, but 

sin^±iasin4r cos f - cos (N + \)a 1 cos (N + \)a 



sin a sin (a/2) 2 sin a sin(a/2) 4 sin (a/2) 2 sin a sin(a/2) 

and the second term oscillates more and more quickly as N increases. So its 
multiplication on any normal function and integration will give a result which 
tends to zero as N -> oo. Therefore l/4sin 2 (a/2) = E^o^n ) 

(cos a) relation 

can be used under integration without any loss of rigor. Now C^(cosa) 
can be replaced by using the addition theorem and we get because of the 
orthonormality property of the spherical functions: 

J sin 2 a/2 n+1 v ' 
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So a 4-dimensional spherical function is indeed a solution of the Eq.41 and the 
corresponding eigenvalue is determined from the po = me 2 / (n + 1) condition 
(m being electron mass), which gives a well known expression for the hydrogen 
atom levels E = — me 4 /2(n + l) 2 , n = 0, 1, . . .. 

Now it should be no surprise that a solution of the Eq.42 is proportional to 
a 5-dimensional spherical function. To prove this, as the previous experience 
suggests, it is necessary to decompose (1 — cos a) -1 into a series of Gegenbauer 
polynomials and afterwards, by using the addition theorem (51), replace them 
by the spherical functions. 

Suppose 1/(1 — x) = E^Lo A]yCffl 2 \x), when because of the orthogonality 
of the Gegenbauer polynomials and of the normalization 



/ 



.2u^(3/2), x l2 , _2(N + 1)(N + 2) 



{l-x 2 )[C$'*\x)Ydx = 



we get 



. 27V + 3 



^- 2WTWT2) / (1 + a;)Cg/2)W ^ 
The integral can be evaluated by using relations 



d 



C$ /2 \x) = —P N+1 (x), P„(l) = 1, and J P n {x)dx = if n ^ 0. 

So 



doc | 



a _ 2N + S 1 _ v> 27V + 3 ( 3/2), , 

^-(ivTWT^j' I^o^"ir (^ + i)(^ + 2) GiV (cosa) " 

8tt 2 

^ f/\r i 1V7U i ^^^™(^)^VnZm(^ ) 5 
Nnlm l 7V + J-A^ + Z ) 



Therefore 



sin 2 a/2 (N+l)(N + 2) 

and eigenvalues of the Eq.42 will be Ajv = (N + 1)(N + 2)m 2 (here m stands 
for the particle mass, not an index of the spherical function), with the corre- 
sponding solution 

$ N (n) = ]T A n i m Y Nn i m (n) . 

nlm 
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Stereographic projection in the Wick-Cutkosky model 



Let us return to Eq.27 and consider again the equal mass case 



m + f p + -P 



m 



+ ^--Pj *(p)=_/dg ^ 



(p - 9) 2 



(52) 



here P = (0, M) - Euclidean 4-vector and M is the bound state mass. 

Equation (52) exhibits explicitly only 50(3) symmetry (because it con- 
tains a fixed 4-vector P), but now we know that a higher symmetry can be 
hidden. To find out this, let us transform the 4-dimensional momentum space 
onto the 5-dimensional sphere surface by the stereographic projection. R.h.s. 
of the equation will be changed at that according to formulas (38) and (40) 

A Ua * (g) - M sec -2^i / ^ (e' 4 /2)1»(g) 

f 2 * (P ~ 1? ~ 16?r 2 2 J sin 2 a/2 ' 

As for the l.h.s. , we will have 



m + f p + -P 



m + f p - -P 



m 4 - -m 2 M 2 + — M 4 + 2m 2 p 2 - -M 2 p 2 + p A + (p ■ P) 2 = 
2 16 2 v 7 

( m 2 -iAf 2 ) 2 + 2 P 2 tg 2 ^l( m 2 -iAf 2 ) + 

p 2 tg 2 ^( P 2 tg 2 ^ + Af 2 cos 2 e 3 ) . 
It is convenient to choose po - 5-sphere radius as: pi = m 2 — (1/4)M 2 , when 



m + Ip + 2^ 



m 2 tg 2 ^c OS 2 e ; 



m + [p — —P 



vl 



2 4 ®4 . 

p sec — + 



2 4 ®4 

= p sec — 



2 M2 • 2^ 2^ 

H sin B4 COS 



and after introduction of the new F(Q) = P(6i, 62, 63, 64) = sec 6 ^^(p) 
function , Eq.52 takes the form 



2 m2 ■ 2 ^ 2 ^ 
p + sin 64 cos 63 



16tt 2 ^ sin 2 a/2 



dn' 



(53) 
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Note that sinG^cosOs = Pa/pq (P is that point, which corresponds to the p 
4- vector in the stereographic projection ), therefore Eq.53 is invariant under 
such rotations of the 5-dimensional coordinate system, which don't change 
the 4-th axis. So Eq.52 possesses a hidden 50(4) symmetry. This symmetry 
was not explicit because ^-momentum space is not orthogonal to the 4-th 
axis (34. But with the inverse stereographic projection, if we choose its pole 
on the axis we get fc-momentum space, which is orthogonal to e±. This is 
schematicly shown in Fig. 10: 




Figure 10: k- momentum space. 



Therefore Eq.52, rewritten in terms of the k- variables, should become explic- 
itly 50(4) invariant. Let us get this new equation. According to the last 
equality in (37) 

■ n n ^ Po~ k2 

sin 64 COS 63 = cos 7 = — = —r -r . 

Po Po + « 



41 



Furthermore, because of Eq.38 

1 



Po 



2 1 2 1_ 

sin 2 a/2 ~ (k - k' f S6C 2 S6C 2 ' 



and according to Eq.40 



dQ' = 4 sec" 8 (- I dk' 
Po \2 



The substitution all of these into the Eq.53, after introduction of the new 
Q(k) = cos ( r f/2)F(Q) unknown function, will give 



But 



• 2 m 2 (pg - k r 

Po 4 (rf + A; 2 ) 2 
4 7 



Po 



sec 



1 + tg 



= 1 + 



and we finally get 



(vl + fc 2 ) 2 + K(pI - fc 2 ) 2 

q Po 



m »pg + ™ (*»)» + 2^1 2 2 
Po V Po 



rn 



Q(k) = 
A /■ Q{k') 



Since the equation just obtained is 50(4) invariant, its solution should have 
the form Q(k) = / (k 2 )Y n i m {Q) , which enables to perform the angular inte- 
gration in Eq.54 and, as a result, to get a one-dimensional integral equation 
for f{k 2 ). Let us denote x = Vk 2 and let a be the angle between k and k' 
4-vectors, then 

1 



(k - k') 2 



x 2 i_ 2 ^cosa+^ 



I x' 2 i_24coscH-4t 



n'=0 



= ^^o(f) n '^ ) (cosa), if x>x' ' 



^^(^"'^(cosa), if z 



(a/) n X 



According to the addition theorem 



it T 1 I'm' 



2tt 2 



+ - 1 - I'm' 
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therefore 



2tt s 



T = £ ~T — 



l\n' 



(k - k 1 ) 



X 



— y6(x - x')+ 



n'+2 



X 



lWm'(^)^ n Vm'(^') 



The substitution of this and = f(x)Y n i m (Q) into Eq.54 will give (because 
of orthonormality of the spherical functions and dk' = l/c'pdl/c'lc^ 4 ^')) 



bo + - 2 ) 2 + ^"- 2 ) 2 



Apt 



2A 
n+ 1 



oo 



,n+3 



^S(x - y) + ^9(y - x) 



x 



y 



Since 



(Po + + ^ - x 2 ) 2 = W + *T 1 1 + '-hi I - = 



2\2 



/to 



,2^2 



4pi 



4p( 



,2\2 



to + ^/ o 

£*0 



PoM W 2 2 p M x 
x\[Pq + x x 



m 



m 



then this integral equation can be rewritten also as 



PoM 



Po + x + x)[p + x ZZ~ X J J\ x ) = 



2Xp 



2 oo 




2/ 



n+3 



m 



9(z - 2/) + ~^®(y - z) 



x n+2 y 



m 2 (n + 1) 

or introducing s = x/po undimensional variable: 

(l + S 2 + M s )(l + s2 _M s) /(s) = 



/to ^ 



2A 



J 



m 2 (n+l) ^ 



r n+3 



B(s - r) + -£r9(r - s) f(r) dr 



(55) 



It is clear from this equation that for the f n (s) = f(s) function, if n ^ 0, we 
have / n (0) = boundary condition at the point s = 0. 

Let us show that Eq.55 is equivalent to a second order differential equation. 
Let 



R(s,r) = 



r 



n+3 



; Q(s - r) + — ^(r - s) 



gn+2 v ' j.n- 
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Since (d/ds)Q(s — r) = 5(s — r), for the derivatives over s we will have: 



d 



r 



n+3 



-R(s, r) = -{n + 2)-^36( S - r) + n^6(r - s) 



n— 1 



^2 j,n+2> 

-^R(s, r) = (n + 2)(n + 3)^e(s - r)+ 



n-2 



ra(n - l)-^ ZT 9(r - s) - 2(n + l)£(s - r) 



Let us try to choose ^4 and 5 coefficients such that 

R(s,r) = -2(n+l)6(s-r) 



d 2 Ad B 
+ - -r + ~z 



ds 2 s ds s^ 
which means the following system 

(n + 2)A - 5 = (n + 2)(n + 3) , + 5 = -n(n - 1) 

with A = 3 and 5 = — n(n + 2) as the solution. Therefore 



d 2 3 d n(n + 2) 



ds 2 s ds 



M 2 ' 

(l + s 2 ) 2 -^s 2 



/(*) = 



4A 7 4A 
J S(s - r)f(r)dr = -— /(s) 



m 



Let us introduce the new variable t = s 2 and the new unknown function 
(f>(t) = t[{l + t) 2 - {M 2 /m 2 )t]f{t). Because 

Id d , d 2 d d d d d d d 2 
- — = 2— and — 7 = —2s— = 2— + 2s— — = 2— + At— . 

s ds dt ds z ds dt dt ds dt dt dt z 



We finally get 

d 2 



dt 2 



Kit) + 



A 



t[m 2 {l + t) 2 - M 2 t] 



n(n + 2) 
At 2 



(56) 



The asymptotic form of this equation (when n ^ 0) in the t — > or t — > oo 
limit is 

d 2 0n _ n(n + 2) _ 
dt 2 4* 2 0n ' 

with solutions £ 1+n / 2 and t~ n l 2 . At the origin n ~ t~ n l 2 behavior is not 
adequate, because then f n ~ t _ ( n / 2 ) _1 ? which contradicts to the / n (0) = 
condition. Let us show that at the infinity, on the contrary, <j) n ~ t 1+n / 2 
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behavior is not good, because then f n ~ £( n / 2 )~ 2 = s n ~ 4 ? which contradicts 
to the integral equation (55). Indeed, when s — > oo, left hand size of the 
equation will be of the order of s n , and the r.h.s. of the order of 



s n+2 



r n+3 f n (r)dr + s ? 



-n— 1 



dr = -s n 2 + 



1 s 

J r n+3 f n (r)dr 



Let sq be some big enough number, so that when r > sq we can use the 
asymptotic form of the solution. Then 

s s s 

j r n+3 f n {r)dr « J r n+3 f n {r)dr + J r 2n ~ l dr = 

s 

i s 

-(s 2 " - 4") + / r" +3 /„(r)dr . 
-1, /r" +3 /»(r)* ~ ,"" 2 . 



Therefore 



So the r.h.s. of the Eq.55 turns out to be of the order of s n ~ 2 and it can't be 
equal to the l.h.s. 

As we see, Eq.56 should be accompanied by the following boundary con- 
ditions 



</>n(t) ~ , if * , 

<f> n (i) ~ r n / 2 , if t oo . 



(57) 



In Eq.56 £- variable changes from to oo. For numeric calculations it 
is more convenient to have a finite interval. So let us introduce the new 
z = (1 — t)/{\ +t) variable, which changes from —1 to 1. Using 

d 2 d d 2 d -2 d 

dt ~ ~(l + t) 2 dz' dtf~ dt(l + t) 2 dz ~~ 



4 



d 



+ 



4 



(1 + t) 3 dz ' (l + £) 4 <k 2 2 l ] dz 4 l ; efe 2 



we get 



{1 ~ Z) ^ + 2{1 - Z) ^z- + 



A 



n(n+ 2) 
m 2 - (M 2 /4)(l-z 2 ) ~ 1-z 2 



0n = , 
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which can be somewhat simplified if we introduce the new g n (z) = (l + z)(l 
z 2 ) n / 2 (j) n {z) function, for which the equation looks like 



[ } dz* 



d 2 9n . dg n 
+ Znz— h 



dz 



A 



m 2_(M 2 /4)(l-z 2 ) 



— n(n + 1) 



9n = 0. (58) 



Eq.57 indicates that at z — > 1 we should have g n ~ (1 — z) n l 2 {l — z) 1+n ^ 2 = 
(l-z) n+ \ and at z -> (l + z) 1+n / 2 (l + z) n / 2 = + Therefore 

the boundary conditions for the Eq.58 are g n (—l) = g n (l) = 0. 



Integral representation method 

Transition from Eq.52 to an one-dimensional integral equation can be per- 
formed also in another way, which is not directly connected with the hidden 
symmetry of the equation. This time let us consider a general case of unequal 
masses and let mi = m + A, = m — A. Eq.27 takes the form 

Hp) = 



TT 2 



((m + A) 2 + (p - z^P) 2 ) ((m - A) 2 + (p + ^ 2 P) 2 )] 1 / dq-^f^ -(59) 



Note, that the solution of the equation 



is, as we have seen earlier 

$(p) - cos 6 ^F^ m (^) - (1 + cos 6 4 ) 3 sin" 6 4 sin 1 e 3 Y lm (G 2 , 6i) 

,2 , _2N-n-3i^ihz //^ ^ \ yimip) 



(^ + m^)-^|p1^ m (9 2 ,e 1 ) = 



(p2 + m 2)n+3 ' 

because 

m 2 -^ 2 2m|p| . |p| 

cos 64 = — ~ ~ , sm 64 = — 7. o ano - sm "3 = l — r 

m z + p z m z + £r |p| 

(Here |p| = is the length of the 4- vector p = (p,Pa) ). 
Using the following parametric representation 

11} 



2_i [(1/2)(A + 5) + (1/2)(A - B)z] z ' 
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we will have 

[(m - A) 2 + (p + i m Pf]- l [{m + A) 2 + (p - ^P) 2 ]" 1 = 



1 r , A 2 A . o o , A. 



m 



(60) 



Here J9 2 , = m2 ~ (1/4)M 2 , as earlier. 

All of this gives a hint to check, how changes, when inserted in the r.h.s. 
of the Eq.59, the following function 



®rdm(p,z) = 



yimip) 



First of all let us evaluate Jdq ^^fi- . According to Eq.29 

J dq ^ 



®nlm{q, Z) 

(p - q) 2 



A 2 A A 

dqyim{q) (1 + —o ~ 2 -z)pI + Q 2 + i(z - —)q ■ P 
m z m m 

l 

[(p " q?Y l = (n + 3)J u n+2 du J dqy lm (q)x 



-(n+3) 



X 



A 



A 



A. 



u{(l + — 2 - 2— + g 2 + z(> )g • P} + (1 - u)(p - q) 



m m 
which can be transformed further as 



-(n+4) 



A 2 A 9 9 / A 
(1 + — - 2-z)p\ + q 2 + z(* - -)g • P 



m 



+ (l-u)(p-q) 2 = 



1 A 

g + )P- (l-u)p 

2 m 



+ u(l - u)p 2 + -u 2 (z - — ) 2 P 2 + 
7 4 m 



A 2 . A 

^2 



wp5(l + — ^ - 2—z) + iu{l - u)(z )p ■ P , 

therefore let us introduce the new k = q+ ^u(z — ^)P — (1 — u)p integration 
variable: 

/ dq ^ m } q q) f = (n + 3) / w n+2 ^ / ^ 3>, m (£ + (1 - u)p)[k 2 + u(l - u)p 2 + 
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\ u \ z - ^)2 P 2 + up 2n + _ 2 ^ z) + m(1 _ _ . p]-(«+4) 

4 m nr m m 

For the solid harmonics the following equation holds 

_^ l k 

yim{a + b)=Y. E M-fc,m-/i(o) * 

4ir(2l + + m)\(l - m)\ 
(2k + 1)(2/ - 2fc + l)(fc + — + m — k — fi)\(l — m — k + n)\ 

Let us decompose 3^m((l ~~ + k) according to this relation and take into 
account that if / ^ 0, then / Yi m (k)d^Q = 0. Therefore only the first term 



^y lm [(i - u)p\y 00 (k) = y lm [(i - u )p\ = (1 - u yy lm (p) 

of this decomposition will give a nonzero contribution. Using also dk = 
(l/2)tdtd^n, where t = k 2 , and J d^Q = 2tt 2 , we will get 



Jdq 



(p - g) 2 



= 7r 2 (n + 3) y/ ro (p) / u n+2 (l - u) l dux 





eft t 



t + u{\- U)p 2 + -M 2 (Z - — ) 2 P 2 + 

7 4 m 



A 2 A A 

+ — ^ - 2—z) + m(l - u)(z )p ■ P 

■■■ m 



nr m 
Over t the integral is of the type 

tdt °r dt 



-(n+4) 



oo 



That is 



n+4 







[* + a] 



n+3 



a 



dt 



[t + a] n + A (n + 2) (n + 3) 



a 



-(n+2) 



/" dq ^ nlm(y(llZ ^ ~ ^ 
{p ~ q) 2 n + 2 



^(p) /(I " «)' 



(1 - n)p 2 + - -) 2 P 2 + 
4 m 



A 2 A A 
p 2 (l + — - 2—z) + z(l - u)(z - —)p ■ P 



-(n+2) 



du 



m m 
or, after introduction of the new t = 1 — u integration variable: 

,2 1 



jdq 



(p — q) 2 n + 2 



yuv) jj\^\p 2 + i{z-^)p-p- 
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4 V m 



A 2 A 1 A x9 2 
— - 2-z) + T (z - -) 2 P 2 



nr m 



-(n+2) 



eft 



Let us denote 



and 



A 2 A A 

a = (1 + — ~ - 2—z)p 2 +p 2 + i(z )p ■ P 

m z m m 



A 2 A , 9 1, A x9 9 
6 = (1 + - 2-z)^ 2 + -(z - -) 2 P 2 



nr m 



m 



We have 

j t l dt 



o 



[t(a-6) + 6]"+ 2 



= (-1) 



t( n + l-Z)! r 

(n + 1)! da'/ [t(a-6) + 6]^+ 2 



1)' 



,n 



+ 1-/)! ^ 



(n + 1)! da' (a -6) 



eft 



* + fc/(a-6)] 



- Mln-/+2 



n 



/)! 5> z 1 [ 1 



+ 1)! da 1 a-b [b n ~ l+1 a n ~ l+1 . 



(-iy 



/)! ^ a 11 - 1 + a n - l - l b + 



a 



n-l-2u2 



b 2 + ... + b 



n—l 



(n + 1)! da 1 



a' 



-l+lfon-l+l 



(-i)' 



/)! ^ 



(n + 1)! da z k=0 



]T a~ {n ~ l+1 ~ k) b~ {k+1) = 



(n-l)\ g (n-fc)! a -(n-k+i) b -(k+i) 



(n + 1)! fc f ( n -Z-fc)! 
where n > / supposition was done. Thus 



jdq 



$nim{q,z) = 2 (n-/)! ^ (n-fc)! 
(p-g) 2 77 (n + 2)! fc ^ (n-/-A;)! X 



A 2 A 



1 



A, 



rf(l + ^-2-z) + ^(z--) 2 P 2 
u m 2 m 4 V m 



-(fc+i) 



x 



but 



yimip) 



A 2 A 



A. 



-(n-fc+l) 



(1 + —2 - 2— 2)Po + P + )P ■ P 



nr m 



m 
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$n-fc-2,ZmO, Z) 



therefore 



{n -I) 



I n-l 



[n-k)\ 



i q {p-qf =% J^T2y-k(n-l-k)\ X 

$n-k-2,lm(P> Z ) 

W + £-^) + \(z-i?P 2 ] k+1 ' 
The result is encouraging, because we have got a linear combination of 
<f> n im functions, and therefore it can be expected that the solution of the 
Eq.59 is expressible through these functions. But, before this conclusion is 
done, we should check that nothing wrong happens after multiplication over 
the 

[(m + A) 2 + - i^PfY^m - A) 2 + (p + i^P) 2 }' 1 ■ 
According to formulas (60) and (29): 

[(m + A) 2 + (p - ir]iP) 2 ]- l [{m - A) 2 + (p + i7] 2 P) 2 ]~ l x 

A 2 A 2 2 ^\ r {n - k+1) 
(1 + — r - 2— z)^ + p 2 + 2(2; )p • P 



m 



m 

1 1 



1 

-(n — A; + l)(n - /c + 2) / eft / drcx 



x(l — x) 



2 

n— fc 



ix(t - z)p • P + (1 + 



-1 
A 2 J, 



A. 

m' 



z(z )p • P — 2 — — z)pq 



A 

I . 

m 



m* m 

-(n-fc+3) 



1 1 



-(n - & + l)(n -k + 2) J dt j dxx 



-1 



x(l — x) ? 



{[1 + £| - 2^(^ + 2 - xz)K + p 2 + z(xt + z - - £)p • P}«-fc+3 



therefore 



[(m + A) 2 + (p - z^P) 2 ]"^™ - A) 2 + (p + ^P) 2 ]" 1 / ^ 

^g(n-Q!(n- fc + 2)!; rft>< 



{p - q) 2 



2 ^{n + 2)\{n-l-k)\J l 
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f x(l - X) n k $ n -k,lm(Pi Xt + Z - Xz)dx 



\k+l 



Note, that when < x < 1, then min(z, t) < xt+z—xz < max(z, t), therefore 
if the z-parameter of the § n i m (p, z) function changes in the — 1 < z < 1 range, 
then — 1 < xt + z — xz < 1 and the above given equation can be rewritten as 

[(m + A) 2 + (p - i^Pfr^m - A) 2 + (p + u^P) 2 ]- 1 / dq ^^ = 

J {p-q) 2 

7T 2 ^ (n - /)!(n - /c + 2)! J J , . . J JA 
— V ) dt dx x(l -x) n k d(x 

^(C — — z + xz) 



So we see that if Q n im{q, z) is substituted in the r.h.s. of the Eq.59, the result 
will be a superposition of the same kind functions. Therefore the solution of 
Eq.59 can be expressed in the form 

n— I J, 

<&nlm{p) = E / 9nli Z )®n-k,lm{p, Z)dz . 
fc=0_! 

Then 

$nim(q,z) 



[(m + A) 2 + (p - zr/ 1 P) 2 ]- 1 [(m - A) 2 + (p + n^P) 2 ] -1 / dg- 



(p- g) 

n— Z 

E 

^=0 



2 



E I dzg^zYim + ^f+ip-i^Pf}-^ 

Km - a) 2 + ( P + ^^p) 2 ]- 1 / = 

r = T(n-z/ + 2)!(n-z/-r-/)! X 
i l i 

/ dzg v nl (z) J dtjdx x(l - x) n -^ r x 

-1 -1 o 

} dC 5((-xt- z + xz)$ n - v - T j m (p, C) 
[rf(l + § - 2^) + I( 2 - A )2p2 ]-' • 



-1 
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Let us denote k = v + r and take into account that 

N N-p N k 

E E = E T, A (",k-v) 

therefore Eq.59 can be rewritten as 

n-l 1 

E 

k=0 



n—l x 

E / g k nl(0$n-k,lm(PX)d( = 



2 gg (n - „ + 2)!(n - f - fcjl / W) / * I *° ^ ~ x) * 



-1 o 



1 ^ $(£ - Xt - Z + Xz)$ n -k,lm(Pi 

- J i [P§(1 + & - 2^) + J(* " ifPif ' 
from this we get the following system of non-homogeneous integral equations 
for the Qni(z) coefficient functions: 



fc _ \ * (n — l — v)\{n - k + 2)\ } j. } , M , n ^ k 

^ = 2 g (n - „ + 2)l(n - / - k)l j * / ^ " ^ 
i 



x 



£(C - art - z + xz)g v nl (z)dz 
>, [pg(l + £ - 2^) + I(, - A) 2 pf-" 

Note that we have interchanged integrations over cfc and g?(\ If now integra- 
tion over dz is further interchanged with integrations over dx and dt, Eq.61 
can be rewritten as (after z <-> C replacement) 

- 2 ( n - i/ + 2)!(n - Z - A;)! X 

/(*,C)<ft(C)dC 



where 



l l 



7(2, £) = J dt J dx x{l - x) n ~ k 5(z -xt-C + xQ = 



-1 

1 



y x(l -x) n k J S(z - xt - C + Ooft 
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Let us introduce instead of t the new y = xt + (1 — x)( integration variable: 

1 x+(l-x)( 

I(z, C) = Jdx{l- x) n - k J 5{z - y)dy = J J (1 - x) n - k 5{z - y)dxdy . 

-x+(l-x)( 

The integration area D = D\ U D2 is shown in Fig. 11. 
It is easy to calculate integrals over D\ and D2. 

1 1 

J J (1 - x) n - k 5(z - y)dxdy = J dy5{z - y) J (1 - x) n - k dx = 

1 C U=£ 

i-C 

1 /• / 1 — y^ n ~ fc+1 



n-k + li- Z - y) \—Q) dV = 



C 

1 /l-Z\ n - fc+1 



G(z - C) 



n - k+ 1 VI - C/ 
and 

C 1 

J J d {1- x) n - k 5{z - y)dxdy = J dy5{z - y) J (1 - xf- k dx = 



i+C 



1 j . . . /l + ^ n - fc+1 



n-Zc+1^ v y Vl + C 
1 /l + z\ n - fc+1 



n - A: + 1 VI + (J 
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Therefore 

'(*>C) = 



1 



n — k + 1 
Let us introduce 



1 - z s 



n-k+l 



G(z - C) + 



1 + 2^ 

.i + C 



n— fc+1 



e(c - *) 



and define the 6-function at the origin as 6(0) = \, then I(z, () = 
and our system of integral equations can be rewritten as 

A * (n — l — v)\(n - k)\ 



n-k+l 



2 ^=0 (n- z/ + 2)!(n- / - k)\ 

{n-k + 2)[R(zX)} n - k+ y nl (Qd( 
i[ ri (l + ^-2A C) + i( C -4)2p 2] *-^' 

In particular, <7^(z) = Pn(^) satisfies a homogeneous integral equation: 



(62) 



0n(*) = 



A 



2(n + l)_^§(l + S-2^C) + |(C-^ 2 



(63) 



Note that R(1,Q = R(-1,Q = 0, if |C| ^ 1. Therefore g n (-l) = g n (l) = 0. 
From Eq.63 A eigenvalues are determined. Inverting A = A(M) dependence, 
we get a mass spectrum M = M(A). 

Let us see, to which second order differential equation is equivalent our 
integral equation (63). We have ^R{z, C) = @(C ~~ z ) I+( ~~ ®( z ~~ Oj^c anc ^ 



-R(zX) = ~j^2^{ z ~ 0- Furthermore 



(1 - , 2 ) 



+ 2zi? 



+ ^ 



dR 

dz 



R 2 = 
R = 



dR 

dz 



R- z 



dR' 

dz , 



i-C : 



,e(z-oe(c-z) 



The r.h.s. differs from zero only at the z = ( point. So 

2 



(1 - , 2 ) 



dRjzX) 
dz 



+ 2zR{zX) 



dR(z, Q 
dz 



-R\zX) = 



l-C 



-e(z-c)e(C-^) 
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Since 



d 



d 



(1 - z 1 )— + 2z(n - 1)— - n(n - 1) i? n (z, C) = 



n(n — 



n-2 



dz 2 
dW 2 



dz 



dz 

+ 2zR^P - R 2 \ + n(l - z 2 )iT~ 



id 2 /? 



and R(z,z) = 1, we get finally 

d 2 



(1 - z 2 )-^ + Mn ~ l)f z ~ n(n - 1) j C) = 

- 2n<5(z - c) - 2 ; e(z - CMC - z) 



(64) 



Using this equality, we can transform Eq.63 into a differential equation (note 
that if /(C) is a normal function, without 5(z — C) type singularities, then 

/Iie(*-c)e(c-*)/(C)<*C = o) 



{l-z 2 )^ + 2nz d ^-n{n+l)g n + 



A 



p 2 ) (l + ^-2^z) + \(z-^) 2 P 2 



g n = o 



(65) 



But 



Po 1 + 



A 5 



A 
• 

m 



I / Ay 

4 V m/ 



2—2! I + - ( Z I P 2 = 



m ^hA ( 1 + 4_ 2 M + i (z _a 

v 4 J \ m z m J 4 m 



rrr 1 + 



A 5 



A 



1 



77T 



2— z --M\l-z l ) 
m ) 4 



therefore Eq.65 can be rewritten also as 



(1 - + 2- 



A 



m\\ + ^-2±z)-\M\\-z 2 ) 



dz 

— n(n + 1) 



9n = 



(66) 



From This equation, when A = 0, we get already known to us Eq.58. It 
turns out that even in a general case of unequal masses Eq.66 can be trans- 
formed into an equation of the Eq.58 type by suitable change of variables. To 
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guess this variable change, it is better to rewrite Eq.66 in the form of Eq.56. 
Introducing (f) n (z) = (1 + z)~ l (l — z 2 )~ n ^ 2 g n (z) function, we get 



(1-^ + 2(1-,) 



d(p r 



dz 

n(n + 2) 

Ml + S " 2£*) " \M\l - z 2 ) ~T^~ 



X 



if n = 



and after the t = j-^ variable change (let us note that = — J- and 



1+z 

(1+t) 4 d 2 i (l+tf d v 



+ 



1 



(ft 2 £ 2 



But 



A 



n(n + 2) 



(i+^) 2 



,A 1-t 2 



-m 
t 



-m 
t 



(1 + t) 2 1 + 



m t 



A' 



M 2 



(f n = 



2-(l-* 2 ) 
m 



H 1 + — +1 + 2t 1 + 



A 5 



my 



m 



-m 



A' 



H 1+ ™ + 1 



A' 



+ 4t- 



A 2 



4^ + m Mi-^ [1 + ((1 + ™ )/(1 -» )M2 



((i + £)/(i- 



A- 



therefore, if one more variable change t = j^f/m £ is made, we get an equation 
of the Eq.56 type: 



d 2 (p r , 
dt 2 



+ 



A/ 1 



a! 

m 2 



n(n + 2) 



£[m 2 (l + If - ((M 2 - 4A 2 )/ (1 - § )) t] 4* 2 



y>„ = 



(67) 

As we see, it is enough to find A = F(M 2 ) spectrum for equal masses. Then 
the spectrum for the unequal masses case can be found from the relation 



A 



1-^ 



= F 



M 2 - 4A T 



(68) 



Note, that t — ► t variable change implies the following transformation for the 
variable z : z 



z— — 



y — 
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Using Eq.64, we can transform Eq.62 also into a system of second order 
differential equations 

(1 - z 2 )^ + 2(n - k)zj- z -( n -k)(n-k + 1) J g k nl (z)+ 

k ( n _/_ z/ )!( n _/ c + 2)! 
h(n-v + 2)\(n-l-k)\ X 
9ni(z) = Q 

[ m 2 (1 + 4!_ 2 A z) _l M 2 (1 _ z 2 )]fc -, + l • 

Variable separation in bipolar coordinates 

There exists one more method for Eq.27 investigation. It is based on the 
following idea: transform integral Wick-Cutkosky equation into a partial 
differential equation and try to separate variables in some special coordinate 
system. 

Transition from the integral to the differential equation can be done by 
using the fact that in a 4-dimensional Euclidean space 



1 



A-^r = -4ir 2 S(x) . (69) 



x 



Let us show that this is indeed correct. Derivatives of \ are singular at the 
x = point. Therefore A^ should be specially defined at this point, for 
example, as 



A 1 1- A 1 

A— T = hm A- 



x 2 e^o x 2 + e 

_ . r> i \ /• 1 

dxi x 2 +e ~~ (x 2 +e) 2 



But 4rz^ = -tJ^ and 



1 d 2 1 8 Sx 2 8e 

+ 



x 2 + e dxidxi x 2 + e (x 2 + e) 2 (x 2 + e) 3 (x 2 + e) 3 
Therefore 

e^O x 2 + e I — oo , if X = 
This suggests that A^ is proportional to —5(x). But 

-86 f . dX „ = - W /\ ,. = -8e. 2 7-^ = -4. 2 



(x 2 + e) 3 ^ (r 2 + e) 3 J (* + e) 
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so Eq.69 holds. 

£)2 £)2 £)2 f)2 

Acting on the both sides of Eq.27 by the operator A p = §^ + §^ + §^ + §^1 
and using \j^2 = ~^ 2 5{p -q),we get 

A p [ml + p 2 + (p 4 - %M) 2 ] [to 2 . + p 2 + (p 4 + ^ 2 M) 2 ]$0) = -4A$(p) , 
or, after the introduction of a new function 

= [™i + P 2 + (P4 - w/iM) 2 ] [to 2 , + p 2 + (p 4 + imM) 2 ]<5>{p) , 

4A 



-#(p) = . (70) 



[m\ + p 2 + (j94 — i?7iM) 2 ] [ml + p 2 + (j>4 + i^M) 2 ] 
So the partial differential equation is found. Now we should care about 
variable separation. It turns out that so called bipolar coordinates can be 
used for this goal. 

On a plane, the bipolar coordinates r, a are defined as follows: a is the 



angle indicated in Fig. 12 and expr = ^. 




Figure 12: Bipolar coordinates. 



Because T\ = \J(x + a) 2 + y 2 and r 2 = \J{x — a) 2 + y 2 , this definition im- 
plies 



x 2 + y 2 - a 2 



1 + a) 2 + ?/ 2 
t = — In — , a = arccos , 

2 (x - a) 2 + y 2 + a) 2 + y 2 ^(x - a) 2 + y 2 

To inverse these relations, let us note that 

i !/ T , — 7-\ x 2 + ?/ 2 + a 2 / 2ax 

cosh r = -( e + e J = , sinh r = V cosh r — 1 = 

2 rir 2 rir 2 



cos a = 



x 2 + y 2 — a 2 
r\T2 
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n 2~ 2a y a u 2flZ 

sin a = v 1 — cos 2 a = and cosh r — cos a = 



r\T2 



r\T2 



So 



a sinh r 



asm a 



2/ 



cosh r — cos a cosh r — cos a 

In a 4-dimensional Euclidean momentum space the bipolar coordinates are 
defined through 



Pi = 



asm a 



coshr — cos a 

a sin a 

1 A = 



sin 6 cos (f , P2 = 



asm a 



cos 6 , P4 = 



coshr — cos a 
a sinh r 



sin 6 sin ip , 



cosh r — cos a cosh r — cos a 

It is convenient to choose the a parameter in such a way to have 

a 2 = ml - j]\M 2 = m\- r&M 2 . 



(71) 



if m = 



mi+m 2 



(as it was assumed so far), this is impossible. But let us 



recall from the BS equation derivation that r\\ and 772 are subject to only one 
771 + 772 condition. In other respects they are arbitrary. So let us demand 
m\ — r] 2 M 2 = m\ — 7/1 M 2 : , which gives 



a = 



M 2 + m 2 — m\ 
' 2M 2 

y/Vjni + m 2 ) 2 - 



M 2 - mj + 777,2 
2M 2 



M 2 ][M 2 - (mi - m 2 y 



2M 



Now Eq.70 should be rewritten in the new coordinates. We have 

9 ^ 9 . ,,. 9 9 a 2 sin 2 a a 2 sinh 2 r 

+ P + (P4 - ^iM) 2 = a 2 + 



+ 



(cosh r — cos a) 2 (cosh r — cos a) 2 



a sinh r 



cosh r — cos a 
a sinh r 



2 , 2 
= a + a 



1 — cos 2 a 



+ a z 



cosh 2 r — 1 



227?iM 



cosh r — cos a cosh r — cos a 
Analogously m\ + p 2 + (734 + if] 2 M) 2 = 
Therefore 



(cosh r — cos a) 2 (cosh r — cos a) 2 
2a 

(a coshr — zr/iMsinhr) . 

-(a coshr + imMsmhr). 

coshr— cos a v IA ' 



[ml + p 2 + (p 4 - M71M) 2 ] [ra 2 . + p 2 + (p 4 + ^ 2 M) 2 ] = 
4a 2 

(a cosh r — t^M sinh r) (a cosh r + 2772 M sinh r) 



6(coshr — cos a) 2 
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Now we should express the Laplacian in the bipolar coordinates. Denoting 

a sin a 
cosh T — COS OL ' 



_ asma ^ 

msh t — ms n ' 



&_ , 1 _ _ 1 / 0^ _ (5^ 



y»2 



Here 5^ is the angular part of the 3-dimensional Laplacian A^ 3 ) = + 
If + 4 Let * = ^ + ir, then | = \ + *») , £ = and 
<|? + ^ = 4^^. So a and r should be expressed through z and z 

1, (j?4 + a) 2 + r 2 1, (z + a){z + a 
r = - m -n = - m 



2 (j94 — a) 2 + r 2 2 (z — a)(z — a) 
As for a, using arccos A = i In [A + v^4 2 — 1], we get 



]?4 + r 2 — a 2 

a = arccos 



yj (pi + «) 2 + r 2 ^j{pA - a) 2 + r 2 

— a 2 z. (z + a)(z — a) 

arccos = = - m — . 

^(z 2 -fl 2 )(z2-a 2 ) 2 (z-a)(z + a) 

It follows from the above given equations that 

da —ia da ia dr a dr a 



dz z 2 — a 2 ' dz z 2 — a 2 ' dz z 2 — a 2 ' dz z 2 — a 2 
Thus 



and 



Therefore 



d da d dr d a ( d . d \ 

dz dz da dzdr z 2 — a 2 \dr da] 

d a ( d . d \ 



i- 



dz z 2 — a 2 \dr da) 
d 2 d 2 4a 2 I d 2 d 2 

+ 7r^T = ^ kt—, 5rhr^ + 



dp 2 ' dr 2 [z 2 — a 2 )(z 2 — a 2 ) \dr 2 ' da 2 1 
(note that + iJ^ = 2 ^^.^ commutes with _ 2 ^ q2 , because z = p^ — ir = 

a sinhr-sinhM = flCoth l£a) R t 
coshr— cosh [ia) 2 / 

(z 2 — a 2 ) (z 2 — a 2 ) = [z + a) (z + a) (z — a) (z — a) = 

4a 4 

[( ft + a) 2 + ^][(p 4 - o) ^ + r^ (coshT _ cosa)2 , 
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and for the Laplacian we finally get 
A (cosh t — cos a) 3 I d 2 ^ d 2 



sin a 



(coshr — cos a)' 



cr sin a 



dr 2 da 2 j cosh r — cos a 



a 2 sin 2 a 



.£(3) 



This could be rewritten in a more convenient form by using 

d 2 ( d 2 n d 1 

—, r sin a = sin a — + 2ctga- 1 

and Eq.45 : 

A = 



(coshr - cos a) 3 / cP_ (4) _ \ 
a 2 \dr 2 ) 



1 



(coshr — cos a) 
So Eq.70 in the bipolar coordinates takes the form 

d 2 ^ A 1 



(coshr — cos a) 



dr 2 
A 



- 6® - 1 



(cosh r — cos a) 

* = 



(a cosh r — ir)\M sinh r) (a cosh r + 2772 M sinh r) 
It is possible to separate variables in this equation. In particular, if we take 

*(r, a, 6, </?) = (cosh r - cos a)/(r)y n / TO (a, 6, </?) . (72) 

and use 5^Y n i m = n(n + 2)Y n i m equation, we get for the function / the 
following equation 



dr 2 



+ 



A 



7 = 0. 

(73) 



(a cosh r — ir\\M sinh r) (a cosh r + ir]2M sinh r) 
To simplify, note that because 

sinh 2r = 2 sinh r cosh r , 1+cosh 2r = 2 cosh r and cosh 2r— 1 = 2 sinh r 
the following holds 

(a cosh r — z?7iM sinh r) (a cosh r + i^M sinh r) = 

2 2 2 2 / \ 

a cosh r + 771772 M sinh r + iMa(r]2 — 771) coshr sinh r = 



-{(a 2 + 771772M 2 ) cosh2r + iaM(rj 2 - 771) sinh2r - (771772M 2 - a 2 )} 



But 



(a 2 + 7 7l 77 2 M 2 ) 2 - % l a l M' l {j] 2 - = 



2„2 11 /f2, 
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a 4 + J]\j]\M^ + a 2 M 2 7]\ + a 2 M 2 T] 2 = (a 2 + r/ 2 M 2 )(a 2 + j]\M 2 ) 
therefore there exists a complex number v such that 



2 2 

m 1 m 2 



cosh z/ = 



1 



17111712 



-(a + T]\T]2M ) and sinhz/ = 



-aM(r]2 - m) 



Besides 



r}\r}<iM 2 - a 2 = f]irj 2 M 2 -m\ + rjlM 2 = j] X M 2 -m\ = 

^(M 2 + m 2 - m 2 2 ) -m\ = ^(M 2 - m 2 - m 2 ) , 
and, using cosh {x + y) = cosh x cosh y + sinh x sinh y, we get 
(a cosh r — 2771 M sinh r) (a cosh r + ir^M sinh r) = 



17lil7l2 



cosh (2r + z/) 



1 



\M 2 -m\ - m 2 ) 



2mim2 

Therefore Eq.73, after the introduction of a new variable a = r + 1, will look 

as 



d(T 2 



+ 



A 



m\m2 cosh cr — t[M 2 — (mi — ztt^) 2 ] 



f = 



(74) 



Its form reveals the already known to us fact that the unequal mass case is 
equivalent to the equal mass problem with m 2 = m\m2 and M' 2 = M 2 — 
(mi — m2) 2 . 

Thus, it is sufficient to consider the following equation (for equal masses 
v = and a = r) 



dr 2 



+ 



A 



m 2 cosh t 



\M 2 



f = 



(75) 



As is clear from the bipolar coordinates definition — oo < r < oo, therefore 
Eq.75 should be supplied with boundary conditions at r — > ±oo. These 
conditions follow from the initial integral equation, but instead we will show 
that Eq.75 is equivalent to the Eq.56. Indeed, let us introduce a new variable 
t = sinh2r + cosh2r = e 2r , then = 4^ + 4£ 2 J^. But (1 + t) 2 = 
(2 cosh 2 r + 2cosh2r sinh2r) 2 = 4 cosh 2 r (cosh 2r + sinh 2r) 2 = 4(cosh 2 r)£, 
so cosh 2 r = ^4r~; anc ^ Eq.75 takes the form 



A 



d 2 f 1 df 

~d£ 2+ ~t~dt + \t[m 2 (l + t) 2 - MH] 



4t 2 



/ = o, 
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which after introduction of a new function = y/if coincides to the Eq.56. 
Therefore Eq.57 implies the following boundary conditions for the Eq.75 

f{r) - e - (n+1)|T| , when r -> ±00 . (76) 

Appearance of the spherical function in Eq.72 indicates that the Wick- 
Cut kosky model possesses the 50(4) hidden symmetry even for unequal 
masses. The variable separation is possible just because of this symmetry. 

Note about references 

Only the sources of this review are presented here. Exhaustive bibliog- 
raphy about the Bethe-Salpeter equation and in particular about the Wick- 
Cutkosky model can be found in [13] • 
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